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SYNOPSIS 

In  this  final  technical  report  we  document,  in  three  chapters,  the  results 
obtained  from  our  AFOSR-sponsored  study  on  a  new  ion  acceleration  using 
beating  electrostatic  waves. 

Unlike  previously  known  methods  of  energizing  plasmas  with  electrostatic 
waves,  and  which  accelerate  only  ions  whose  initial  velocities  are  above  a  certain 
threshold  (close  to  the  waves  velocity),  the  new  mechanism  can  accelerate  ions 
with  arbitrarily  small  initial  velocities  (as  illustrated  in  Fig.  (1)).  This  results 
in  significant  improvements  to  the  acceleration/heating  efficiency  (by  as  much 
as  45%)  over  previous  techniques.  The  benefits  to  the  state  of  the  art  of  space 
propulsion  stem  from  the  high  efficiency  of  the  new  mechanism,  its  thersholdless 
nature  and  its  capabilities  of  producing  high  ion  energies  and  its  electrodeless 
character. 


Figure  1:  Calculation  from  first  principles  (using  Hamiltonian  formulation)  of 
ion  acceleration  by  two  beating  electrostatic  waves  (red  curve)  compared  to  the 
effect  of  a  single  wave  (blue  curve)  .  Unlike  the  single  wave  case,  the  2-wave 
interaction  accelerates  ions  whose  initial  velocity  is  below  the  classical  thresh¬ 
old  (dotted  horizontal  line)  to  significantly  high  velocities.  The  thresholdless 
nature  and  efficiency  of  this  acceleration  mechanism  make  it  most  promising  for 
propulsion  applications. 

In  summary  we  have  accomplished  the  following: 

1.  Demonstrated  theoretically  that  a  new  highly  efficient  ion  acceleration 
mechanism  by  beating  electrostatic  (ES)  waves  is  real  and  is  a  fundamen¬ 
tal  phenomenon  that  occurs  only  under  specific  conditions  that  were  not 
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known  before.  We  derived,  from  first  principles,  the  fundamental  neces¬ 
sary  and  sufficient  conditions  for  the  acceleration  mechanism  to  occur. 

2.  We  studied  using  Monte  Carlo  simulations  the  effects  of  collisions  (that 
would  occur  in  a  real  low  temperature  plasma)  on  the  ion  acceleration 
mechanism  and  found  that  collisions  can  enhance  the  effects  due  to  scat¬ 
tering  of  ions  out  of  the  forbidden  region  of  phase  space. 

3.  We  designed  and  constrducted  a  dedicated  experiment  to  study  the  new 
ion  acceleratio  mechanism.  We  were  able  to  launch  and  beat  ion  cyclotron 
waves  in  an  argon  plasma  produced  by  a  helicon  source  and,  using  a  ded¬ 
icated  state-of-the-art  LIF  diagnostic  system  (acquired  through  a  recent 
DURIP  grant),  we  were  able  to  provide  the  first  laboratory  de- 
mostration  of  the  existence  of  this  new  ion  acceleration  mecha¬ 
nism. 

This  final  technical  report  consists  of  three  chapters  summarized  below. 

In  Chapter  1  the  conditions  under  which  a  magnetized  ion  can  be  acceler¬ 
ated  through  a  nonlinear  interaction  with  a  pair  of  beating  electrostatic  waves 
are  explored.  The  electric  field  of  the  beating  waves  can,  under  some  condi¬ 
tions,  accelerate  ions  from  arbitrarily  low  initial  velocity  in  stark  contrast  with 
the  well-known  nonlinear  threshold  criteria  for  ion  acceleration  by  a  single  wave. 
A  numerical  investigation  of  the  fundamental  dynamics  led  to  the  identification 
of  critical  points  on  the  Poincare  cross-section.  A  second-order  perturbation 
analysis  was  carried  out  to  approximate  the  location  of  these  critical  points  and 
define  the  domains  of  allowed  and  forbidden  acceleration.  It  is  shown  that  for 
an  ion  to  be  significantly  energized,  the  Hamiltonian  must  be  outside  the  en¬ 
ergy  barrier  defined  by  the  location  of  the  elliptic  and  hyperbolic  critical  points. 
Despite  the  restriction  on  the  Hamiltonian,  an  ion  with  arbitrarily  low  initial 
velocity  may  benefit  from  this  acceleration  mechanism.  The  resulting  domain 
of  allowed  acceleration  is  significantly  larger  than  that  of  ion  acceleration  by 
either  single  or  non-beating  waves. 

In  Chapter  2  a  numerical  model  of  the  nonlinear  interaction  between  beat¬ 
ing  electrostatic  waves  and  magnetized  ions,  including  collisions,  is  presented. 
Previous  studies  of  the  beating  electrostatic  waves  (BEW)  interacting  with  a 
single  ion  showed  the  ability  of  this  mechanism  to  accelerate  ions  from  arbitrar¬ 
ily  low  initial  velocities,  and  have  revealed  the  fundamental  conditions  for  this 
interaction  to  occur.  The  present  study  extends  the  analysis  to  a  large  number 
of  ions  and  includes  ion-ion  collisions.  The  numerical  investigation  combines  a 
dynamical  description  for  the  ion-wave  interaction  and  a  Monte  Carlo  simulation 
of  the  collisions.  Despite  the  thermalization  role  of  collisions  BEW  acceleration 
was  found  to  yield  larger  heating  rate  and  higher  particle  energies  than  the 
better  known  interaction  with  the  single  electrostatic  wave  (SEW). 

In  Chapter  3  we  report  on  our  experiments  that  demonstrate  the  exis¬ 
tence  of  the  new  acceleration  mechanism.  We  first  report  on  the  excitation  and 
propagation  of  Electrostatic  Ion  Cyclotron  (EIC)  waves  in  an  rf-sustained  argon 
plasma  and  our  measurements  of  the  sought  dispersion  relations.  Such  waves 
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were  used  to  study  the  beating  waves  ion  acceleration  mechanism  as  reported  in 
the  last  section  of  the  chapter.  The  waves  are  excited  by  an  antenna  consisting 
of  two  parallel  metal  plates  inserted  at  the  edge  of  a  plasma  column  with  their 
surface  normal  perpendicular  to  magnetic  field.  The  plates  can  be  driven  in  or 
out  of  phase.  The  in-phase  configuration  couples  better  to  plasma.  It  is  shown 
that  EIC  waves  launched  at  frequencies  between  uci  and  10wcj  propagate  with 
little  damping  at  an  angle  between  82°  and  86°  with  respect  to  the  magnetic 
field.  The  amplitude  of  the  excited  waves  can  be  optimized  by  properly  match¬ 
ing  the  impedance  of  the  driving  circuit  to  the  plasma  and  choosing  the  right 
plasma  conditions.  The  dispersion  relation  was  measured  using  a  phase-delay 
technique  and  was  found  to  be  in  good  agreement  with  the  theoretical  EIC  dis¬ 
persion  relation  over  a  wide  range  of  frequencies.  We  conclude  that  chapter  and 
the  report  with  the  first  laboratory  demonstration  of  the  existence  of 
the  new  ion  acceleration  mechanism  by  directly  measuring  a  significant 
(up  to  40%)  increase  in  the  perpendicular  ion  energy. 
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grant. 


Chapter  1 


Ion  Acceleration  by  Beating 
Electrostatic  Waves: 

Domain  of  Allowed 
Acceleration 


1.1  Introduction 

Stochastic  heating  of  a  magnetized  ion  by  a  single  propagating  electrostatic  (ES) 
wave  has  been  extensively  studied  [1,  2,  3,  4].  Using  first-order  perturbation 
theory  Karney  [3,  5]  was  able  to  derive  analytical  expressions  approximating 
overall  nonlinear  dynamics  of  an  ion  interacting  with  a  single  ES  wave.  That 
work  revealed  the  existence  of  a  threshold  for  the  initial  ion  energy  below  which 
the  particle  cannot  gain  net  energy  from  the  ES  wave.  Skiff  et  al.  validated  these 
findings  experimentally  [6].  The  threshold  is  essentially  the  lower  bound  of  a 
nonlinearly  broadened  resonance  condition  between  the  ion  velocity  and  wave 
phase  velocity.  Ions  with  initial  velocities  below  this  threshold  will  not  gain  net 
energy  and  their  motion  will  remain  coherent,  while  ions  moving  slightly  faster 
than  the  wave  will  be  accelerated  stochastically.  Thus  the  threshold  separates 
two  regions  of  phase  space:  a  Tegular  (or  coherent)  motion  region  of  low  energies 
below  the  threshold  and  a  stochastic  one  -  above  the  threshold.  Nonlinear  ion 
acceleration  by  a  single  wave  is  therefore  always  a  stochastic  process. 

In  1998  Benisti  et  al.  described  a  new  mechanism  for  nonlinear  ion  acceler¬ 
ation  by  ES  waves  [7,  8] .  The  scheme  requires  pairs  of  ES  waves  that  obey  a 
beating  criterion:  their  frequencies  must  differ  by  an  integer  number  of  the  ion 
cyclotron  frequency.  In  essence  the  beating  of  two  or  more  ES  waves  creates  a 
propagating  electric  field  structure  having  a  low-frequency  amplitude  envelope. 
Slow  ions  are  first  accelerated  coherently  through  interaction  with  this  envelope 
and  are  then  further  accelerated  stochastically  by  the  higher  harmonics  of  the 
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beat  wave.  Under  such  conditions  the  single- wave  theory  threshold  disappears 
and  regular  and  stochastic  regions  of  phase  space  become  connected  allowing 
ions  with  arbitrary  small  initial  velocities  to  obtain  high  energies  through  co¬ 
herent  acceleration  followed  by  stochastic  energization.  Consequently,  this  non¬ 
linear  interaction  may  result  in  a  more  efficient  acceleration  mechanism  than  is 
possible  from  the  interaction  with  a  single  wave.  This  new  acceleration  scheme 
has  been  advanced  as  a  possible  explanation  for  the  ionospheric  ion  acceleration 
observed  during  the  Topaz  3  [9]  rocket  experiments.  An  acceleration  mechanism 
that  promises  to  energize  a  larger  portion  of  the  ion  distribution  function  may 
be  promising  to  many  applications  where  the  efficiency  of  ion  heating  is  of  prime 
importance,  such  as  in  spacecraft  plasma  propulsion. 

A  preliminary  numerical  exploration  [10]  of  this  mechanism  revealed  that 
there  are  many  cases  for  which  the  beating  criterion  does  not  lead  to  ion  ac¬ 
celeration.  This  hinted  to  the  possibility  that  the  criterion  is  necessary  but  not 
sufficient. 

In  this  study  we  attempt,  using  analytical  and  numerical  techniques,  to 
define  the  necessary  and  sufficient  conditions  for  the  acceleration  of  the  magne¬ 
tized  ion  through  nonlinear  interaction  with  a  pair  of  propagating  electrostatic 
waves.  In  particular,  we  use  second-order  perturbation  theory  and  numerical 
integration  to  analyze  the  nonlinear  dynamics  by  finding  the  location  of  the 
critical  points  [11]  of  the  motion  on  a  Poincare  cross-section  [12].  These  critical 
points  were  not  observed  in  previous  studies  [7,  8]  since  they  were  limited  to  the 
analysis  of  a  single  trajectory  (i.e.  single  initial  condition)  in  a  given  Poincare 
cross-section.  The  location  of  these  critical  points  allows  us  to  define  the  initial 
conditions  for  the  Hamiltonian  and  the  ion’s  velocity  required  for  a  magnetized 
ion  to  be  energized  by  a  pair  of  beating  waves. 

In  section  1.2  we  overview  the  analytical  formulation  of  the  problem.  In 
order  to  give  context  to  our  study  we  review  in  section  1.3  the  better-known 
dynamics  of  a  magnetized  ion  interacting  with  a  single  ES  wave.  This  is  followed 
in  section  1.4  by  an  overview  of  the  fundamental  features  of  the  new  mechanism 
that  relies  on  beating  waves.  In  the  remaining  sections  we  seek  expressions  for 
the  location  of  the  critical  points  of  the  motion  which  amount  to  a  definition  of 
the  domains  of  allowed  and  forbidden  acceleration. 


1.2  Analytical  Formulation 

We  start  by  defining  the  coordinate  axis  as  shown  in  Fig.  1.1.  Here  we  have 
a  single  ion  of  mass  m  and  charge  q  in  a  constant  and  homogeneous  magnetic 
field,  Bz.  This  ion  interacts  with  a  packet  of  electrostatic  waves  that  propagate 
in  the  positive  x  direction.  Since  we  take  the  waves  to  be  purely  electrostatic 
the  wave-number  fc,  is  parallel  to  the  electric  field  £)  of  each  of  these  waves. 
The  dynamics  of  the  single  ion  in  Fig.  1.1  is  described  by  the  following  equation 


FINAL  REPORT  9/1/04:  ION  ACCEL.  WITH  BEATING  WAVES 


8 


Figure  1.1:  A  single  ion  of  charge  q  and  mass  m  in  a  constant  homogeneous 
magnetic  field  Bi  interacts  with  a  spectrum  of  electrostatic  waves  whose  wave- 
number  and  electric  field  direction  is  parallel  to  the  x-axis. 


of  motion  [13,  7]: 

^  £<  sin(k Xx  -  ujit  +  pi),  (1.1) 

i 

qB/m  is  the  ion  cyclotron  frequency  and  pt  is  the  phase  angle  of 
The  corresponding  Hamiltonian  for  the  system  is  [7] 

H  =  p2 /2  +  —  cos(«;jpsin  6  —  uxr  +  pi).  (1.2) 

In  writing  Eq.  (1.2)  we  have  used  the  fact  that  the  system  is  periodic,  and 
transformed  the  Hamiltonian  into  normalized  action-angle  coordinate  system 
[14],  where  k,  =  fcj/fci,  v,  =  oJi/ujc,  r  =  u>ct,  ex  =  (k  1qEi)/{mul),  p2  =  X2  +  X2, 
and  X  =  kjx,  X  =  dX/dr,  so  that  X  =  psin 9,  X  =  pcosQ.  The  action-angle 
coordinate  system  is  a  special  case  of  polar  coordinates  [12].  In  our  context  0 
corresponds  to  the  cyclotron  phase  angle  measured  clockwise  from  the  y-axis, 
as  indicated  on  Fig.  1.1,  while  p  is  the  normalized  Larmor  radius  which,  in  a 
constant  magnetic  field,  represents  the  normalized  velocity  (perpendicular  to 
the  magnetic  field)  of  the  magnetized  particle  undergoing  cyclotron  motion  in 
the  xy  plane. 

Benisti  et  al.  [7]  defined  a  criterion  for  particle  acceleration  by  multiple  ES 
waves.  They  showed  that  for  regular  and  stochastic  regions  to  be  connected  it 
is  necessary  (but,  as  we  shall  see,  not  sufficient)  to  have  at  least  one  pair  of  ES 


where  u>c  = 
each  wave. 
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waves  such  that 


u>2  ~  =  nuic,  (1.3) 

where  n  is  an  integer.  They  also  report  that  acceleration  is  more  vigorous  for 
n  <  2,  therefore  for  the  sake  of  simplicity,  we  limit  our  analysis  to  the  case  of  a 
single  pair  of  beating  waves,  such  that  n  =  1.  In  addition,  Ref.  [7]  reports  that 
the  maximum  acceleration  is  achieved  when  all  waves  are  of  the  same  amplitude, 
£j  =  £j  =  e.  We  also  set  /c*  =  Kj  —  sq  to  simplify  our  analysis,  and  since  the 
phase  angles,  (pi,  do  not  play  a  fundamental  role  in  this  acceleration  process  [7] 
we  set  all  < pi  =  0.  With  these  simplifications  the  Hamiltonian  (1.2)  becomes 

H  =  p2  [2  +  e[cos(psin0  —  U{t)  +  (1.4) 

cos(psin0  —  i/j-r)]. 

This  Hamiltonian  represents  two  coupled  oscillators:  one  is  the  gyrating  ion  and 
the  other  corresponds  to  the  beating  ES  waves.  We,  therefore,  can  interpret  e 
as  a  coupling  parameter  between  the  two  oscillators. 

Two  approaches  are  taken  in  this  study  for  analyzing  the  system  above. 
First,  we  derive  the  analytical  solution  to  equation  (1.4)  by  applying  a  second- 
order  perturbation  technique  in  conjunction  with  Deprit  modified  [15,  16]  Lie 
transformations  [12,  13],  and  then  compare  it  with  the  results  of  numerical 
integration  of  equation  (1.4).  To  obtain  the  numerical  solution  we  use  the  fourth- 
order  symplectic  integration  algorithm  (SIA4)  derived  by  Candy  and  Rozmus 
[17]  who  showed  its  superiority  over  fourth  order  Runge-Kutta  algorithm  for 
Hamiltonian  periodic  problems,  such  as  ours. 

A  convenient  way  of  representing  both  numerical  and  analytical  solutions  is 
to  plot  the  resulting  trajectories  on  a  Poincare  cross-section  [12].  To  construct 
a  Poincae  cross-section  from  the  numerical  integration  of  Eq.  (1.4)  we  plot  the 
point  intersections  of  the  ion  trajectory  in  three  dimensions  ( p ,  0,  r)  with  the  p-0 
plane  at  specific  time  intervals.  For  integer  values  of  v  this  reduces  to  plotting  p 
vs.  0  at  r  =  2nj,  where  j  —  0, 1, 2, ...  is  a  non-negative  integer.  For  non-integer 
values  of  v  precaution  must  be  taken  for  proper  accounting  of  intersection  points. 
Since  the  magnetic  field  is  constant,  the  normalized  cyclotron  radius  p  is  a  direct 
measure  of  the  perpendicular  ion  velocity.  Therefore  Poincare  cross-sections  give 
direct  visual  insight  into  the  acceleration  process. 

The  visual  interpretation  of  Poincare  cross-sections  is  straightforward.  Ran¬ 
dom  point  distribution  corresponds  to  stochastic  motion  while  regular  patterns, 
such  as  lines  and  ellipses,  will  tell  us  that  the  ion  dynamics  is  analytical  (or 
regular).  For  example,  if  the  wave  amplitude,  e,  is  zero,  Eqs.  (1.1)  and  (1.4) 
reduce  to  a  simple  harmonic  oscillator  and  for  irrational  values  of  v  its  Poincare 
cross-section  shows  a  set  of  horizontal  lines,  indicating  constant  velocity  (which 
corresponds  to  a  free  ion  gyrating  in  a  constant  magnetic  field).  Each  of  these 
lines  represents  an  invariant  of  motion  for  a  given  set  of  initial  conditions  [12]. 
When  the  coupling  parameter  e  is  not  zero,  we  can  treat  the  ion  motion  as  a 
perturbation  of  these  invariants. 
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The  detailed  derivation  of  the  analytical  solution  for  a  particle  interacting 
with  a  single  wave  can  be  found  in  Ref. [5].  However,  a  more  generalized  solution 
for  multiple  waves  is  obtained  [12,  13,  18]  through  Deprit’s  modified  Lie  trans¬ 
formation  in  Refs.  [7,  8].  The  resulting  autonomous  Hamiltonian  derived  from 
Eq.  (1.4)  for  a  non-integer  value  of  v  to  the  second  order  in  the  perturbation, 
e,  is 


where, 


H  =  e{  JVi  ( p )  cos  (ui9)  +  JVi  (p)  cos  {vjO)} 
+  e2{S?(p)  +  S?(p) 

+  Sei','i(p)cos[(uj-Ui)e}}, 


S?(p)  = 


°6 


(P)  = 


+ 


J_  y'  rnJm(p)Jmip) 

2  p  ^  Vi—m  ’ 

'  m——oc 

J_/  A  rnJm{p)Jl._Vi+m{p)  ^ 

2  p  ^  Vi-m 

Nm=— oo 


E 


m=  —  oo 


m  Jm  (  P)  J Vi  -  Vj 
Vj-m 


(1.5) 


(1.6) 


Jm  is  the  Bessel  function  of  the  first  kind  of  order  m,  and  J'  represents  the 
derivative  of  the  Bessel  function  with  respect  to  its  argument.  When  v,  is  an 
integer,  the  summations  are  performed  over  all  m  /  Ui  to  avoid  singularities. 
When  v  /  integer,  the  first  order  terms  in  Eq.  (1.5)  disappear  [7],  and  the 
equation  becomes  more  tractable. 

The  Hamiltonian  in  Eq.  (1.5)  is  autonomous  and  therefore  itself  is  an  in¬ 
variant  of  motion.  Curves  of  constant  H  in  a  Poincare  cross-section  represent 
the  complete  analytical  solution  of  the  problem  to  second  order.  We  now  wish 
to  compare  the  Poincare  cross-sections  of  the  analytical  solution  to  those  ob¬ 
tained  through  numerical  integration  of  Eq.  (1.4).  It  should  be  noted  that  we 
should  not  expect  to  see  any  stochastic  behavior  on  the  Poincare  cross-sections 
obtained  from  the  analytical  solution. 

As  with  most  phase  diagrams,  critical  points  define  the  dynamics  of  motion. 
Since  the  system  is  not  dissipative  we  expect  to  find  two  types  of  critical  points: 
elliptic  and  hyperbolic.  As  we  shall  show  later,  the  location  of  the  critical 
points  is  the  key  to  determining  which  initial  conditions  lead  to  acceleration  or 
trapping.  The  task  before  us  is  to  find  these  critical  points. 

We  now  explore  particle  dynamics  as  a  function  of  wave  amplitude  and 
frequency,  and  in  terms  of  the  location  of  critical  points  on  the  Poincare  cross- 
section.  Using  numerical  solutions  we  will  demonstrate  that  when  critical  points 
are  absent  in  the  regular  region  of  the  Poincare  cross-section  (as  in  the  single 
wave-particle  interaction),  the  particle  will  not  gain  net  energy. 
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1.3  Single  Wave  Interaction 


0  400  800  1200  1600 

T 

Figure  1.2:  Time  evolution  of  the  normalized  velocity  p  for  a  particle  interacting 
with  a  single  off-resonance  wave  (e  =  10,  u  =  24.3).  The  threshold  derived  in 
Ref.[5]  and  given  by  Eq.  (1.7)  represents  the  boundary  between  the  regular  and 
stochastic  domains  and  is  shown  as  a  horizontal  dashed  line.  This  picture  is 
qualitatively  similar  to  the  case  of  the  non-beating  waves. 

In  this  section,  in  order  to  create  a  context  for  our  study,  we  summarize 
the  results  obtained  by  Karney  [5,  3]  for  the  interaction  of  a  magnetized  par¬ 
ticle  with  a  single  electrostatic  wave.  When  the  wave  frequency  is  exactly  an 
integer  number  of  the  ion  cyclotron  frequency  it  is  said  to  be  an  on-resonance 
wave,  otherwise  it  is  an  off-resonance  wave.  In  both  cases  Karney  found  that  a 
threshold  given  by 


p=v--Je,  (1-7) 

separates  the  regions  of  regular  and  stochastic  motion  [5,  19],  as  shown  Fig.  1.2, 
which  illustrates  for  the  numerical  solution  of  the  equation  of  the  motion  (1.1)  a 
typical  case  (e  =  10,  v  =  24.3).  Below  the  threshold,  indicated  by  the  horizontal 
dashed  line  ( p  ~  20),  we  observe  that  the  ion  motion  is  regular  and  we  can 
predict  its  behavior  well  by  means  of  perturbation  theory.  More  importantly, 
as  long  as  the  ion’s  initial  velocity  is  in  that  region  it  is  clear  that  the  ion 
will  not  gain  net  energy  from  the  wave  [5].  When  the  initial  velocity  is  above 
the  threshold,  the  ion  moves  stochastically  and  eventually  gains  net  energy,  as 
shown  by  the  upper  trajectory  in  Fig.  1.2.  Therefore,  in  the  case  of  interaction 
with  a  single  wave  the  ion  gains  energy  only  chaotically  when  its  initial  velocity 
po  exceeds  the  threshold. 
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Figure  1.3:  Time  evolution  of  the  normalized  velocity  p  for  a  particle  interacting 
with  two  beating  waves  (e  =  10,  i>i  —  24.3,  Vj  =  25.3).  Two  particle  trajectories 
are  shown,  one  with  an  ion  trapped  below  p  <  v  —  ^Je  threshold  and  another 
resulting  in  an  accelerating  ion.  For  both  trajectories  po  are  the  same  but  0O 
are  different. 

This  is  an  important  point  and  we  will  contrast  it  later  with  the  case  of 
beating  waves  where,  as  we  shall  see,  the  two  regions  of  phase  space  can  become 
connected  under  some  conditions,  and  a  particle  with  an  arbitrarily  small  initial 
velocity  in  the  regular  region  can  be  accelerated  through  the  threshold  to  high 
values  of  p. 

1.4  Multiple  Wave  Interaction 

One  of  the  early  investigations  of  two  ES  waves  interacting  with  ions  was  done 
by  Chia  et  al.  [13,  20]  who  conducted  an  analytical  study  of  the  interaction 
to  first  order  in  the  perturbation  and  found  that  for  the  waves  with  even-even 
or  odd-odd  combinations  of  frequencies  v,  there  existed  a  vertical  separatrix 
allowing  “infinite”  heating.  In  1998  Benisti  et  al.  [7,  8,  9]  showed  that  to 
observe  particle  acceleration  through  the  threshold  boundary  one  needs  to  use 
perturbation  theory  to  at  least  second  order.  According  to  their  analysis,  the 
normalized  condition  for  the  new  acceleration  mechanism  already  given  by  Eq. 
(1.3)  is 

Vi-Vj=n,  (1.8) 

where  n  is  an  integer.  This  amounts  to  a  wave  beating  condition.  By  analyzing 
a  single  ion  trajectory  obtained  from  numerical  solutions  of  Eq.  (1.1)  Benisti  et 
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al.  showed  that  as  long  as  condition  (1.8)  is  satisfied  the  ion  will  gain  significant 
energy  from  the  waves. 

Subsequently,  we  performed  numerical  investigation  [10]  based  on  the  same 
single  trajectory  method  and  found  that  some  ion  trajectories  did  not  lead  to 
ion  acceleration  even  if  condition  (1.8)  was  satisfied,  as  shown  in  Fig.  1.3.  We 
concluded  that  to  find  the  necessary  and  sufficient  conditions  for  ion  acceleration 
we  need  to  examine  the  complete  Poincare  cross-section  and  find  the  critical 
points  of  motion,  as  will  be  done  in  sections  1.5  and  1.6. 

We  first  investigated  the  case  of  two  electrostatic  waves  when  condition  (1.8) 
does  not  hold.  The  picture  is  qualitatively  very  similar  to  that  of  the  single  wave- 
particle  interaction  shown  in  Fig.  1.2.  Ions  with  energies  below  some  threshold 
maintain  coherent  motion  and  do  not  gain  net  energy.  This  picture  can  change 
drastically  if  condition  (1.8)  is  satisfied. 

As  we  already  mentioned  above,  the  case  of  the  off-resonance  beating  waves 
(both  u\  and  iq  are  not  integers)  provides  us  with  much  more  tractable  analytical 
expressions.  Consequently,  we  will  focus  our  attention  on  such  cases.  Some 
qualitative  analysis  of  the  on-resonance  beating  waves  is  still  possible  and  will 
be  done  in  section  1.7. 

In  Fig.  1.4  we  show  typical  Poincare  cross-sections  obtained  by  numerical 
integration  for  iq  =  24.3  and  iq  =  25.3.  The  panels  in  this  figure  illustrate 
the  effect  of  increasing  wave  amplitude.  The  phase  diagram  consists  of  two 
regions,  stochastic  and  regular,  just  as  for  the  single  wave  interaction.  However, 
unlike  the  single  wave-ion  interaction,  the  two  regions  are  “connected”.  By 
“connected”  we  mean  that  an  ion  with  low  initial  velocity  can  undergo  first 
regular  and  then  stochastic  acceleration,  reaching  high  energies. 

For  low  perturbation  strength  (low  values  of  e)  the  regular  region  extends  to 
values  of  p  approximately  predicted  by  Eq.  (1.7).  However,  as  e  is  increased  the 
regular  region  quickly  shrinks  to  the  vicinity  of  the  elliptic  critical  point  (desig¬ 
nated  E  on  Fig.  1.4.)  Notice  that  the  elliptic  point  is  located  at  pe  ~  v/2  and 
Be  =  Jr.  In  dimensional  quantities  the  position  translates  into  relation  between 
the  ion  velocity  and  the  wave  phase  velocity,  vioTl  ~  —0.5u/k.  Eventually,  as 
the  wave  amplitude  is  raised  above  values  shown  on  Fig.  1.4,  chaotic  motion 
dominates  the  phase  diagram. 

We  now  gauge  how  well  the  second-order  perturbation  analysis  compares  to 
the  numerical  solutions.  Fig.  1.5  indicates  a  good  degree  of  agreement  between 
the  two.  Even  though  the  detailed  structure  of  the  regular  motion  lines  is  not 
captured  with  the  analytical  solution,  the  latter  does  predict  the  position  of  the 
lower  elliptic  (E)  as  well  sis  the  hyperbolic  point  (H)  rather  well.  On  the  other 
hand,  our  analytical  approach  breaks  down  in  the  stochastic  region,  as  should 
be  expected.  Therefore  the  critical  points  shown  by  the  analytical  solution  to 
be  at  p  >  25  in  Fig.  1.5  (which  can  be  said  to  describe  a  “homoclinic  tangle”  or 
“stochastic  layer”)  are  in  reality  covered  by  the  stochastic  motion,  as  shown  by 
the  numerical  solution. 

However,  as  described  in  Ref.  [8],  even  in  that  region  of  phase  space  the 
overall  ion  motion  could  be  approximated  by  first-order  orbits,  for  small  e. 
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Figure  1.4:  Poincare  cross-sections  showing  numerical  solutions  for  a  particle 
interacting  with  two  beating  off-resonance  waves  (vz  =  24.3,  Vj  =  25.3).  The 
stochastic  region  occupies  a  greater  fraction  of  the  phase  space  as  the  wave 
amplitude  is  increased. 
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Figure  1.5:  Poincare  cross-section  for  a  particle  interacting  with  two  beating  off- 
resonance  waves  (e  —  10,  Vi  =  24.3,  Vj  =  25.3).  a)  Analytical  solution  showing 
the  existence  of  hyperbolic  and  elliptic  points  marked  by  H  and  E  respectively, 
b)  Numerical  solution  also  showing  the  locations  of  the  critical  points. 
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Figure  1.6:  Analytical  Hamiltonian  as  a  function  of  p  at  9  =  7r.  The  elliptic 
an  the  hyperbolic  points  corresponding  to  the  ones  shown  in  Fig.  1.5  axe  the 
minimum  and  the  maximum  of  H.  (vi  —  24.3,  t/j  =  25.3). 


1.5  Topology  of  the  Phase  Diagram 

As  with  any  phase  diagram,  each  curve  on  the  Poincare  cross-section  corresponds 
to  a  given  set  of  initial  conditions.  In  the  case  of  a  particle  interacting  with 
beating  waves  we  are  mainly  concerned  with  the  hyperbolic  and  elliptic  critical 
points  designated  H  and  E  respectively  on  Fig.  1.5.  It  is  clear  by  tracing 
trajectories  in  Figs.  1.4  and  1.5  that  an  ion  with  Hamiltonian  lying  between 
the  Hamiltonian  values  corresponding  to  points  E  and  H  does  not  gain  net 
energy  from  the  waves  i.e.  does  not  reach  the  stochastic  region  where  it  can 
be  vigorously  accelerated.  Instead  the  corresponding  phase  space  trajectories 
circulate  around  the  elliptic  critical  point  E  or  cover  the  full  range  of  cyclotron 
phase  angles  (0  <  9  <  7r)  while  remaining  below  H. 

Other  features  of  the  Poincare  cross-section  in  Fig.  1.5  worth  mentioning 
(but  not  central  to  our  discussion  on  acceleration)  are  the  two  degenerate  saddle 
points  that  can  be  seen  at  the  intersection  of  the  fourth  curve  (outward)  from 
the  elliptic  point  with  the  0-axis  and  the  “primary”  (upper)  separatrix  that  is 
approximately  near  the  eleventh  curve  (outward)  from  the  elliptic  point. 

It  is  relevant  to  note  in  this  context  that  the  Hamiltonian  of  various  tra¬ 
jectories  increases  monotonically  from  the  Hamiltonian  value  at  point  E  to  its 
value  at  point  H,  as  shown  in  Fig.  1.6.  The  figure  shows  the  Hamiltonian  as 
a  function  of  p  for  0  =  7r  and  illustrates  that  the  location  of  the  elliptic  and 
the  hyperbolic  points  could  be  found  by  determining  the  local  minimum  and  the 
maximum  of  H. 
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Therefore,  for  given  values  of  i >i  and  e,  the  inequality 

He  <  H(p0, 90)  <  Hh,  with  pQ  <v  -  y/l  (1.9) 

defines  the  forbidden  acceleration  domain,  where  H e  and  Hh  are  the  Hamilto¬ 
nian  values  for  the  elliptic  (E)  and  the  hyperbolic  (H)  points,  and  subscript  0 
refers  to  initial  conditions.  By  “forbidden  acceleration  domain”  we  mean  here 
the  domain  of  initial  conditions  for  which  an  ion  cannot  reach  the  stochastic 
region  of  phase  space  where  it  can  be  vigorously  energized.  All  other  ion  tra¬ 
jectories  then  lie  in  the  allowed  acceleration  domain  of  phase  space.  The  ions 
in  the  allowed  acceleration  domain  will  be  affected  by  the  waves  strongly.  The 
restriction  on  po  in  Eq.  (1.9)  is  needed  because  ions  with  po  >  v  —  s/t  will  not 
be  trapped  in  the  energy  barrier  between  the  elliptic  and  the  hyperbolic  points 
(i.e.  in  the  forbidden  acceleration  domain),  as  shown  in  Fig.  1.6. 

The  “trapping”  criterion  in  Eq.  (1.9)  given  in  terms  of  the  Hamiltonian 
should  be  contrasted  with  the  threshold  criterion  for  interaction  with  a  single 
wave  given  by  Eq.  (1.7).  It  is  clear  that  unlike  the  single- wave  case,  an  ion  with 
initial  velocity  p0  below  the  “threshold”  can  still  be  accelerated  to  high  energies 
if  the  corresponding  Hamiltonian  is  outside  the  range  described  by  Eq.  (1.9) 

From  the  point  of  view  of  plasma  acceleration  one  would  like  to  limit  the 
number  of  particles  trapped  in  the  forbidden  acceleration  domain  [He  <  H[p0,  60)  < 
Hh )■  The  rest  of  the  ions  gain  much  higher  energies  through  first  regular  (if 
their  initial  energy  is  low)  and  then  stochastic  acceleration,  as  shown  in  Figs. 

1.4  and  1.5.  However,  even  the  trapped  particles  can  escape  into  the  stochastic 
domain  if  we  consider  a  collisional  plasma,  given  that  a  trapped  particle  may 
gain  enough  energy  during  a  collision  to  overcome  the  energy  barrier,  Fig.  1.6, 
as  we  have  shown  through  particle  simulations  reported  in  Ref.[21]. 


1.6  Critical  Points 

To  define  the  domains  of  allowed  and  forbidden  acceleration  described  by  Eq.  (1 .9) 
we  need  to  find  the  location  of  the  critical  points  E  and  H.  We  now  seek  ana¬ 
lytical  expressions  for  both. 

Since  both  points  are  the  extrema  of  the  Hamiltonian,  that  task  can  be 
achieved  by  setting  the  time  derivative  of  p  and  9  to  zero  simultaneously  [12, 11]. 
Utilizing  Hamilton’s  equations  of  motion  in  conjunction  with  Eqs.  (1.5)  and  (1.6) 
we  get 


p  = 

dH 

—  =  e{uiJVi(p)  sm[ui6) 

(1.10) 

+ 

VjJvAp)  sin(^6>)} 

+ 

e2(i/«  -  V3)Se"Vi(p)  sin[(^  -  i>i)0\  =  0, 

e  = 

=  £{J'v.[p)  cosine) 

(1.11) 

+ 

J'u.[p)  cos [vjd)} 
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+  e2{S'^(p)  +  S^(p) 

+  S™(p)  cos[(j/i  -  Vj)d}}  =  0. 

When  both  wave  frequencies  are  off-resonance  {vi,Vj  /  integer),  the  equations 
above  simplify  because  the  first-order  terms  drop  out,  and  we  are  able  to  obtain 
the  position  of  critical  points  analytically. 

For  V  /  integer,  the  S/(p)  term  in  Eq.  (1.6)  could  be  simplified  to  an 
algebraic  equation  containing  only  few  Bessel  functions  [13]: 

-  Jn-M J-(„-,)W].  (1.12) 

As  a  result  of  this  simplification  we  can  reduce  the  Sg"1'3  (p)  term  down  to 

S?” ( P )  =  f  ST' (P)  +  77s i  (P)-  (113) 

I'i  Vj 

The  details  of  this  derivation  are  given  in  the  appendix.  We  can  therefore 
express  the  Hamiltonian  (1.5)  in  terms  of  the  simplified  S/  (p)  function  only, 

H  =  e2 1  (1  +  —  cos[i2,  -  l/j]0)Si‘  (p) 
l  vi 

+  (1  +  —  cos[i/<  -  i'j]0)Sii  (p)|-  (1.14) 

vj  > 

As  we  will  show  later  in  this  section  our  analysis  breaks  down  for  small  values 
of  v.  Consequently  we  take  iij>  1,  and  remembering  that  n  —  1  with  p/vj  = 
p/{vi  + 1)  ~  p/ u/1  —  1/vf  + ...)  we  can  approximate  p/i/,  ~  p/vj.  Dropping  the 
subscripts  in  S’/  (p)  we  have 


1  +  —  cos(j/j  —  Vj)9 

Vi 


sVi(p)  +  s^(p) 


(1.15) 


Finally,  we  substitute  Eq.  (1.12)  for  each  of  the  S/  (p)  functions.  Expressing 
everything  in  terms  of  we  get 

H  =  -4^—  (l  +  -cose)  (1.16) 

8  sin  U7T  \  v  ) 

—  jfi/_i(p)  j_(1/_i)(p)  +  jl/(p)j-u(p) 

+  Jv+\{p)J-(v+\/p)  —  Jv+2{p)  "J-(i/+2)(p)l  1 


where  we  have  replaced  i/t  with  u. 
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We  are  now  ready  to  find  the  position  of  points  E  and  H.  Prom  Eq.  (1.10) 
as  well  as  from  Figs.  1.4  and  1.5  we  see  that  these  points  lie  at  0  =  7r.  This 
reduces  equations  (1.10)  and  (1.11)  to 

^-(1  -  -)L(p)  =  0,  where  (1.17) 

op  V 

•^'(Pi*/)  =  —  Ju—  (r— l)(p)  d"  Jv(p)J—v{p) 

+  Jv+\(p)J-{v+l){p)  -  Jv+2{p)J ~(v+2){p)- 

Prom  Eq.  (1.17)  we  can  express  p  as  a  function  of  Lip)  and  L'(p),  where  the 
prime  denotes  the  derivative  with  respect  to  p,  and  arrange  the  resulting  ex¬ 
pression  as 


P 

v 


F(P,  v) 


1 

1 


1  L{p,v) 

v  L'(p,  v)  ’ 

£  1  l(p^) 

v  v  L’(p,  v) 


or 


=  0. 


(1.18) 

(1.19) 


The  first  and  second  roots  of  Eq.  (1.19),  pE  and  pn,  for  a  given  value  of  u, 
correspond  to  the  locations  of  the  elliptic  and  hyperbolic  points  respectively. 

The  solution  to  Eq.  (1.18)  for  a  range  of  v  is  shown  in  Fig.  1.7,  while  the 
F{p,v )  of  Eq.  (1.19)  is  plotted  as  a  function  of  p/v  for  v  e  [55.001,55.999]  in 
Fig.  1.8.  It  is  now  important  to  discuss  the  behavior  of  these  solutions. 

As  we  discussed  above,  we  are  only  considering  the  cases  with  u  /  inte¬ 
ger.  Unfortunately,  due  to  the  asymptotic  behavior  of  the  Bessel  function  near 
the  “turning  point”  [22],  defined  as  ptp  =  \/v(y  +  1),  the  solution  approaches 
different  limits  as  v  gets  close  to  an  integer  from  different  sides.  This  peculiar 
behavior  is  demonstrated  for  the  case  of  v  ~  55  in  Fig.  1 .8.  Therefore,  no  simple 
analytical  expression  could  be  obtained  for  the  position  of  the  hyperbolic  point. 
Instead,  Fig.  1.7  shows  a  range  of  solutions  for  the  location  of  this  point.  On  the 
other  hand,  the  elliptic  point  is  very  well  defined.  It  is  seen  that  for  sufficiently 
large  values  of  v  the  location  of  the  elliptic  point,  according  to  our  second-order 
perturbation  analysis  is  at 


Pe  —  2 1  Qe  =  v.  (1-20) 

Although  no  simple  expression  could  be  found  for  the  location  of  the  hyperbolic 
point,  we  see  from  Fig.  1.7  that  its  location  asymptotes  (at  large  values  of  v)  to 
a  value  of  p/v  ~  (0.8  —  1.0),  therefore  we  may  approximate 

pH  ~  0.9i^,  9h  =  7 r.  (1.21) 

It  is  also  important  to  mention  that  because  of  the  asymptotic  behavior  of 
Bessel  functions,  the  elliptic  and  hyperbolic  points  could  not  always  be  found 
for  small  values  of  v.  For  these  cases  while  the  expression  for  the  Hamiltonian, 
Eq.  (1.16),  is  still  valid,  our  analysis  of  the  allowed  and  forbidden  acceleration 
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Figure  1.7:  Solution  to  Eq.  (1.18)  as  a  function  of  v.  The  elliptic  and  hyperbolic 
points,  E  and  H,  could  not  always  be  found  for  small  u.  For  v  >  9  both  points 
could  be  found  for  any  values  of  v.  That  region  is  marked  as  “continuous”  on 
the  figure. 


Figure  1.8:  Plot  of  the  function  F(p,v),  given  by  Eq.  (1.19),  showing  how  the 
second  root  (corresponding  to  the  hyperbolic  point)  moves  away  from  p/v  =  1  as 
v  changes  from  55.001  to  55.999.  The  first  root  on  the  left,  marked  E,  represents 
the  elliptic  point  and  the  subsequent  roots  represent  the  hyperbolic  points  for 
the  various  values  of  pjv. 
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domains  does  not  apply.  The  range  of  v  for  which  we  could  always  find  both 
critical  points  is  designated  as  “continuous”  in  Fig.  1.7. 

Another  limitation  is  placed  on  our  analysis  by  its  independence  on  e  in 
locating  the  critical  points,  as  seen  from  Eq.  (1-17).  Extensive  numerical  explo¬ 
ration  of  the  weak  dependence  of  the  critical  point  locations  on  e  suggests  the 
following  corrections  to  expressions  (1.20)  and  (1.21),  which  are  reminiscent  of 
the  e  dependence  in  the  single  wave  interaction  [5],  Eq.  (1.7), 

PE  -  - — 2^’  ^B  =  7r>  (i-22) 

pH  ^  v-y/e,  6h-tt.  (1.23) 

The  physical  interpretation  and  significance  of  these  two  points  could  be 
understood  as  follows.  Looking  back  at  Fig.  1.1  we  see  that  both  points  cor¬ 
respond  to  the  particles  moving  180°  out  of  phase  with  the  electric  field  of  the 
wave.  We  can  note  from  Fig.  1.6  that  at  the  hyperbolic  point  the  ion  velocity  is 
vh  —  —ui/k  +  vtr  and  the  Hamiltonian  is  at  the  maximum,  while  at  the  elliptic 
point  the  ion  velocity  is  ve  =  (—u>/k  +  vtr)/2,  and  the  Hamiltonian  is  at  the 
minimum.  The  difference  in  the  Hamiltonian  of  the  two  points  then  forms  an 
“energy  barrier”,  which  an  ion  must  overcome  to  be  accelerated  by  the  wave. 

Another  way  to  understand  the  importance  of  the  elliptic  point  is  to  realize 
that  at  this  point  the  energy  exchange  between  the  ion  and  the  waves  is  minimum 
and  the  situation  is  equivalent  to  stable  equilibrium  for  a  pendulum.  Any  small 
perturbation  from  that  equilibrium  will  only  cause  small  oscillations  about  it. 
This  implies  that  in  the  immediate  neighborhood  of  point  E  the  ion  energy 
cannot  be  altered  sufficiently  to  push  the  ion  into  the  stochastic  region,  as 
seen  in  Fig.  1.5.  On  the  other  hand,  the  hyperbolic  point  corresponds  to  the 
unstable  equilibrium  of  the  pendulum  and  any  small  perturbation  from  it  will 
cause  significant  changes  in  the  ion  motion,  i.e.  escape  into  the  stochastic  region 
and  subsequent  vigorous  heating. 

1.7  Beating  Waves  (On-Resonance) 

When  we  choose  vt  and  Vj  to  be  both  on-resonance,  the  overall  behavior  becomes 
much  more  complicated  and  no  simple  analytical  expression,  as  in  the  previous 
section,  can  be  found.  Fig.  1.9  shows  a  typical  example,  the  case  with  e  =  10, 
Vi  =  24,  and  vj  =  25.  One  of  the  major  differences  with  respect  to  the  off- 
resonance  case  is  that  now  we  have  two  hyperbolic  points  which  do  not  lie  at 
6  =  vr.  Nevertheless,  we  could  still  find  their  positions  by  solving  equations 
(1.10)  and  (1.11)  numerically. 

The  location  of  the  elliptic  point  is  easier  to  obtain  since  the  first-order  terms 
in  Equations  (1.10)  and  (1.11)  drop  out  for  9  =  7r,  as  in  the  off-resonance  case. 
Also,  Fig.  1.9  shows  that  the  analytical  solution  exhibits  much  more  complicated 
chains  of  the  critical  points  at  large  values  of  p.  While  in  reality  these  points 
are  “inside”  the  stochastic  region,  this  graphical  picture  illustrates  why  the 
analytical  treatment  of  the  on-resonance  case  is  more  challenging.  However,  we 
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note  that  the  locations  of  both  the  elliptic  and  the  hyperbolic  points  even  in 
this  case  are  very  close  to  those  determined  by  Eqs.  (1.6).  Indeed,  by  studying 
numerically  the  solutions  to  Eq.  (1.4)  for  both  on-resonance  and  off- resonance 
cases  we  conclude  that  the  locations  of  elliptic  and  hyperbolic  points  for  both 
cases  could  be  well  predicted  by  Eqs.  (1.6). 

1.8  Summary  and  Concluding  Remarks 

The  beating  criterion  (u>i  —  u>2  =  tujc )  proposed  by  Benisti  et  al.  [7,  8]  can 
allow  a  magnetized  ion  to  be  energized  by  a  pair  of  beating  electrostatic  waves. 
The  importance  of  this  mechanism  stems  from  its  ability  to  accelerate  ions  with 
arbitrarily  low  initial  velocity.  It  has  become  clear  however,  (see  Fig.  1.3)  that 
this  criterion  is  not  sufficient  for  acceleration. 

In  order  to  better  define  the  criteria  for  acceleration  we  investigated  multiple 
ion  trajectories  (multiple  initial  conditions)  on  the  same  Poincare  cross-section. 
This  analysis  led  to  the  identification  of  critical  points  on  the  phase  diagram. 
Vigorous  ion  acceleration  now  can  be  explained  in  terms  of  the  location  of  these 
points  in  the  region  of  regular  motion.  A  second-order  perturbation  analysis  of 
the  equation  of  motion  allowed  us  to  derive  the  criterion  defining  the  allowed 
and  forbidden  acceleration  domains  in  terms  of  the  location  of  these  points. 

According  to  this  analysis,  for  a  pair  of  beating  {v\  —  U2  —  1),  electrostatic 
waves  interacting  nonlinearly  with  a  magnetized  ion,  significant  ion  acceleration 
can  occur  as  long  as  the  Hamiltonian  of  the  system  does  not  satisfy  the  following 
“trapping”  criterion 

He  <  H{p0;90)  <  Hh,  with  p0  <  i>  -  y/e,  (1.24) 

(which  strictly  applies  when  v  1).  If  the  ion’s  initial  conditions  do  not 
satisfy  the  above  trapping  criteria,  the  ion  can  be  accelerated  from  arbitrarily 
low  initial  velocity  through  the  region  of  regular  motion  to  the  stochastic  region 
where  substantial  energization  can  occur. 

Regular  ion  acceleration  is  a  much  slower  process  than  stochastic  energiza¬ 
tion  [10].  However,  as  the  wave  amplitude  is  increased,  the  region  of  stochastic 
motion  can  extend  down  to  low  initial  velocities.  It  is  important  to  note  that 
the  trapping  criterion  is  in  terms  of  the  (initial)  Hamiltonian  and  not  just  the 
(initial)  velocity  ( po ). 

The  necessary  (wave-beating)  condition  stated  in  Eq.  (1.8)  along  with  the 
avoidance  of  the  trapping  criterion  stated  in  Eq.  (1.24)  represent  two  necessary 
and  sufficient  conditions  for  the  beating-wave  ion  acceleration  mechanism  to 
occur. 

Finally,  it  is  important  to  mention  that  the  trapping  criterion’s  independence 
of  the  wave  energy  e  is  a  consequence  of  the  second-order  nature  of  the  analysis. 
In  light  of  the  albeit  weak  dependence  on  e  in  the  single-wave  criterion  in  Eq. 
(1.7),  we  investigated  numerically  over  a  wide  range  of  non-dimensional  param¬ 
eters  (e  «  5  —  100,  v  «  10  —  50)  location  of  the  critical  points  and  found  that 


Figure  1.9:  Poincare  cross-section  for  a  particle  interacting  with  two  beating 
on-resonance  waves  showing  a  more  complicated  picture  than  that  of  the  off- 
resonance  case  shown  in  Fig.  1.5.  (e  =  10,  =  24,  u3  =  25).  a-1)  and  a-2) 

analytical  solution,  b)  numerical  solution. 
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the  same  dependence  on  applies.  Therefore,  given  the  locations  of  the  critical 
points  in  Eq.  (1.6)  we  can  rewrite  the  trapping  criterion  in  Eq.  (1-24)  as 

H[( v  -  y/i)/2\  7 r]  <  H(p0;  90)  <  H{v  —  vr), 
with  pa  <  v  —  y/i. 

While  the  above  study  offers  insight  into  the  fundamental  problem  of  a  sin¬ 
gle  ion  interacting  with  two  beating  waves,  the  relevance  of  the  mechanism  to 
practical  problems  involving  a  plasma  rests  on  resolving  a  number  of  issues:  1) 
the  effects  of  oblique  wave  propagation,  as  recently  studied  in  Ref.  [23];  2)  the 
effects  of  wave  dispersion;  3)  the  extension  to  collection  of  particles  and  the 
role  of  collisions.  This  last  effect  was  a  subject  of  numerical  investigation  in 
Ref.  [21].  In  that  work  we  found  that  collisional  scattering  enhances  ion  ener¬ 
gization  by  providing  an  escape  mechanism  for  the  ions  trapped  in  the  forbidden 
acceleration  domain  of  phase  space. 


1.9  APPENDIX:  S6U  j(p)  Term  Simplification 


Using  Eq.  (1.8)  and  substituting  for  Jm(p)  and  J'm{p)  with  the  following  iden¬ 
tities  [22] 


9rn 

Jjn-l(p)  +  Jm+l(p)  ~  - Jm(p), 

P 


Jm-l{p)  Jm+l(p)  —  Jm(P)’ 


we  can  rewrite  Sg"1'3  (p)  as 


gUi '"j  (p)  —  i  ^ Jm+2(p)Jm+l{p) 


+  Jm(p)Jm+l{p)  _  y' 

u  —  m.  4-—* 


-  E 
+  E 

-  E 

Now  we  use  identity  [24] 


v  —  m 
Jm{p)  Jm— 1  (p) 

V  —  m 

«An+ 1  2  (p) 

V  —  771 

Jm  (p)  Jm+1  (p)  1 
v  —  m  j 


+  E 


-E 


V  —  771 

Jm— 1  (p)*An-f2(p) 
v  —  m 

Jm—\  2(p) 

v  —  m 

Jm{p)Jm— 1  (p) 


E  Jm+p(p)Jm(p)  =  — 

u  —  m  sin  7T 


v  —  m  sm  TTi/ 

m=  —  oo 

which  is  valid  for  p  >  0  to  simplify  Eq.  (A-26)  to 


srj(p) 


[2J„+2(p) •/_(„+!)  {p) 


8sin7ri/ 


(A-25) 


(A-26) 
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which,  with  the  help  of  identities  (A-25),  may  be  easily  shown  to  equal  to 


v  sm  irv 


S?*(p )  = 

~Ji/—l{p)J—(v—l)(p)] 

+  7 — ~~tt  — - 7 - — r  [Jv(p)J—i/(p) 

(l/  +  l)sm7r(l/+l)L 

~Jv+2(p)J-(v+2)(p)\- 


Chia  et  al.  [13]  showed  that  Sp  (p)  can  be  simplified  as 

—  Jvi  - 1  (p)  J-  (Vi  - 1)  (p)]  • 

It  is  then  clear  that 

S?"  ( p )  =  £■%*  (p)  +  (P).  (A-27) 

Vi  Vj 

Finally,  we  caution  that  the  relation  (A-27)  holds  only  for  the  special  case 
of  L>i  /  integer  and  Uj  —  i/*  +  1. 


Chapter  2 

Effects  of  Ion  Collisions  on 
Ion  Acceleration 


2.1  Introduction 

Acceleration  of  magnetized  ions  by  beating  electrostatic  waves  (BEW)  is  a  non¬ 
linear  phenomenon  that  may  be  occurring  in  nature  and  may  have  interesting 
applications  to  various  problems  including  spacecraft  propulsion.  Observations 
made  with  the  Topaz  3  rocket  [9]  indicated  that  ions  are  accelerated,  in  a  region 
of  natural  electrostatic  wave  activity,  in  the  topside  ionosphere  to  the  escape 
velocity.  A  puzzling  issue  is  that  initial  velocities  of  these  ions  are  significantly 
below  the  previously  known  threshold  required  for  resonant  acceleration  by  elec¬ 
trostatic  waves.  The  threshold  was  derived  in  the  context  of  ion  interaction  with 
a  single  electrostatic  wave  (SEW)  [3,  5]. 

Benisti  et  al.  [7,  8]  proposed  a  non-resonant  acceleration  mechanism  that 
relies  on  nonlinear  interaction  of  an  ion  with  a  pair  of  beating  waves.  They 
showed  that  if  the  criterion 


nuc  =  w  2  — k>i,  (A-l) 

is  satisfied  between  any  pair  of  electrostatic  waves,  ions  can  be  accelerated 
from  an  arbitrary  low  initial  velocity.  Equation  (A-l)  states  that  the  difference 
between  the  frequencies  of  the  two  beating  waves  ai\  and  u>2  should  equal  to  an 
integer  multiple,  n,  of  the  ion  cyclotron  frequency  toc.  This  is  in  great  contrast 
with  the  well  known  SEW-ion  interaction  studied  theoretically  by  Karney  et  al. 
[3,  5],  Zaslavsky  et  al.  [24,  4],  and  Chia  et  al.  [13,  20],  and  experimentally  by 
Skiff  et  al.  [6].  These  studies  showed  that  for  acceleration  to  take  place  ion 
initial  velocity  has  to  be  within  a  resonance  band  of  the  wave  velocity.  Another 
fundamental  understanding  obtained  from  these  studies  was  that  the  ion  motion 
during  SEW-ion  acceleration  is  always  stochastic. 

Choueiri  and  Spektor  [10]  investigated  the  beating  wave  acceleration  mech¬ 
anism  theoretically  and  found  that  while  Eq.  (A-l)  is  necessary  for  the  accel- 


26 


FINAL  REPORT  9/1/04:  ION  ACCEL.  WITH  BEATING  WAVES 


27 


eration  to  occur,  it  is  not  a  sufficient  condition.  Spektor  and  Choueiri  [25] 
derived  and  verified  the  necessary  and  sufficient  conditions  for  interaction,  also 
discussed  in  section  2.2.  When  these  conditions  are  satisfied,  an  ion  with  an 
arbitrary  low  initial  velocity  can  accelerate  through  a  regular  (non-stochastic) 
motion  in  the  electric  field  of  the  beating  waves,  then  reach  a  threshold  above 
which  acceleration  continues  more  vigorously  (stochastically). 

Since  the  BEW  acceleration  can  increase  the  perpendicular  velocity  of  all 
ions,  as  opposed  to  only  the  resonant  part  of  the  distribution  function,  it  is  of 
particular  interest  to  spacecraft  propulsion  applications  where  acceleration  or 
heating  efficiency  is  of  prime  importance.  In  order  to  obtain  the  first  indication 
of  the  existence  of  this  mechanism  we  have  designed  and  built  a  dedicated 
experiment  using  a  helicon  source  and  RF  antenna  to  launch  pairs  of  beating 
waves  [26]. 

In  order  to  guide  the  design  of  the  experiment  and  help  in  interpreting  its 
results  we  needed  a  model  of  the  interaction  that  more  resembles  the  case  of  a 
real  plasma  than  does  the  single  ion  model.  Such  a  model  should  account  for  the 
interaction  of  the  waves  with  large  amount  of  particles,  and  most  importantly 
ion-ion  collisions.  In  this  study  we  present  such  a  model  based  on  using  Monte 
Carlo  techniques  to  describe  collisions,  and  solving  the  equation  of  motion  be¬ 
tween  collisions.  We  use  the  simulation  to  study  parametrically  the  effects  of 
ion  collisions  on  the  heating  rate  and  attainable  average  energy  for  both  SEW 
and  BEW. 

In  section  2.2  we  review  the  collisionless  model  that  describes  the  interac¬ 
tion  of  a  single  particle  with  a  spectrum  of  electrostatic  waves.  We  also  review 
previous  findings  resulting  from  that  model.  In  section  2.3  we  present  the  nu¬ 
merical  model  that  allows  tracking  a  large  number  of  ions  and  account  for  a 
finite  collision  rate.  In  section  2.4  we  present  and  discuss  the  results  of  our  nu¬ 
merical  investigation,  and  in  section  3.7  we  summarize  our  findings  and  deduce 
a  phenomenological  picture  that  illustrates  the  fundamental  differences  between 
BEW  and  SEW  ion  acceleration. 


2.2  Single  Particle  Model 

A  theoretical  model  for  beating  electrostatic  waves  interacting  with  a  single 
ion  is  given  in  [7,  8,  25].  The  description,  which  latter  in  this  chapter  is  aug¬ 
mented  with  inclusion  of  collisions  and  the  ability  of  tracking  many  ions,  is 
shown  schematically  in  Fig.  2.2.  The  schematic  shows  an  ion  in  a  constant 
magnetic  field  Bi  and  electrostatic  wave  traveling  in  transverse  direction  x. 
The  wave  interacts  with  a  gyrating  ion  causing  a  change  in  its  Larmor  radius. 
Because  the  magnetic  field  is  constant,  an  increase  in  the  Larmor  radius  di¬ 
rectly  corresponds  to  the  increase  in  the  ion’s  perpendicular  velocity  and  thus 
its  kinetic  energy.  The  equation  of  motion  governing  the  interaction  between 
a  spectrum  of  propagating  electrostatic  waves  and  a  single  ion  can  be  easily 
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► 


Figure  2.1:  A  single  ion  of  charge  q  and  mass  m  in  a  constant  homogeneous 
magnetic  field  Bi  interacts  with  an  electrostatic  wave.  The  wavenumber  and 
electric  field  of  the  wave  is  parallel  to  the  z-axis. 


derived  [3,  13]: 

d?x 

^  E*  s>n(ki*  -  Wit  +  (fii),  (A-2) 

i 

where  x  and  t  are  the  coordinate  and  the  time  variables,  q  and  m  are  the  charge 
and  the  mass  of  the  ion,  wc  =  qB/m,  Ei,  ki,  wl,  and  (pi  are  the  amplitude,  wave 
number,  frequency  and  phase  of  the  iih  electrostatic  wave.  While  Fig.  2.2  shows 
a  single  wave,  a  similar  picture  can  be  drawn  for  a  spectrum  of  electrostatic 
waves  traveling  in  the  same  direction.  It  is  convenient  to  normalize  the  above 
equation  and  express  it  in  the  canonical  form  [14,  12]: 

H  =  p2/ 2  +  ^  —  cos(Kipsin  9  —  ViT  +  (fii) ,  (A-3) 

l 

where  H  is  the  Hamiltonian  of  the  system,  =  kj/kj ,  Ui  =  Wi/wc,  t  =  u>ct,  Ei  = 
( kiqEi) / (mw 2),  p 2  =  X2  +  X2,  and  X  —  kix,  X  =  dX/dr,  so  that  X  = 
psinO,  X  =  pcosO.  Where  9  is  the  cyclotron  rotation  angle  measured  clockwise 
from  the  y-axis  as  shown  in  Fig.  2.1  ,  and  p  is  the  normalized  Larmor  radius. 
Equations  (A-2)  and  (A-3)  could  be  solved  numerically  with  either  conventional 
4th  order  Runge-Kutta  scheme  or  a  symplectic  approach.  We  have  used  the 
symplectic  integration  method  developed  by  Candy  and  Rozmus  [17]  to  study 
the  behavior  of  a  single  ion  interacting  with  one  or  two  propagating  electrostatic 
waves. 

We  were  able  to  confirm  [25]  that  while  a  single  electrostatic  wave  produces 
some  ion  heating  under  restricted  (resonance)  conditions,  two  beating  waves 
can  result  in  ion  acceleration  from  arbitrary  low  initial  velocities.  We  have 
also  shown  that  Eq.  (A-l)  describes  the  necessary,  but  not  sufficient  condition 
for  that  heating  to  take  place.  The  necessary  and  sufficient  conditions  for  ion 
heating  by  the  beating  electrostatic  waves  are  [25]: 

Tl(jJc  —  U)2  —  ^1 1 

H{p;  6)  =  Hh  >  H(p~  v-  i/e\9  =  7r), 


(A-4) 

(A-5) 
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where  Hh  is  the  value  of  the  Hamiltonian  evaluated  at  the  hyperbolic  point  as 
described  in  Ref.[25]  and  shown  schematically  in  Fig.  2.2. 

The  schematic  shows  the  possible  acceleration  processes.  Particle  1  is  ac¬ 
celerated  stochastically  in  both  cases.  Particle  2  with  initial  energy  below  the 
SEW  resonance  threshold  (p  =  u  —  s/e)  is  affected  by  BEW  but  is  never  al¬ 
lowed  to  reach  the  stochastic  region.  While  particle  3  remains  unaffected  by 
the  SEW  interaction,  it  can  be  effectively  accelerated  by  BEW  through  regular 
(non-stochastic)  motion  that  allows  it  to  reach  stochastic  region,  where  more 
rigorous  acceleration  takes  place. 


Figure  2.2:  Poincare  cross-section  (phase  diagram)  schematic.  The  figure  shows 
typical  trajectories  for  various  initial  conditions  of  ion  interacting  with:  a)  a 
single  electrostatic  wave  (SEW),  b)  beating  electrostatic  waves  (BEW). 
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2.3  Including  Collisions 

Collisions  alter  the  picture  described  in  the  previous  section  drastically.  Without 
collisions  an  ion  whose  initial  velocity  (Hamiltonian)  is  below  that  corresponding 
to  Eq.  (A-5),  particle  2  in  figure  2.2b  will  never  be  effectively  accelerated  by  the 
waves.  However,  a  collision  would  instantaneously  change  that  ion’s  trajectory 
and  place  it  in  a  part  of  phase  place  where  Eq.  (A-5)  is  satisfied. 

In  this  section  we  consider  ion  collisions  only.  Coulomb  ion  collisions  are  of 
interest  since  they  thermalize  the  heavy  species  energy  and  since  our  main  focus 
is  ion  heating. 

To  introduce  collisions  into  our  numerical  model  we  follow  the  classical  work 
of  Takizuka  and  Abe  [27] .  We  model  Coulomb  collisions  as  small  angle  binary 
collisions  and  assume  that  on  a  sufficiently  small  time  scale  we  can  uncouple 
particle  motion  from  collisions.  Thus  our  algorithm  consists  of  two  parts.  We 
move  all  particles  between  collisions  according  to  the  equation  of  motion  pre¬ 
scribed  by  the  single  particle  collisionless  model,  Eq.  (A-2).  We  then  use  the 
Monte  Carlo  approach  to  determine  randomly  the  collision  partners  and  the 
scattering  angles  for  each  collision. 

2.3.1  Overall  implementation 

1.  We  first  choose  a  time  step  At  smaller  than  the  ion-ion  relaxation  time 
calculated  at  the  initial  temperature  of  the  ions. 

2.  Using  a  4th  order  Runge-Kutta  scheme  we  then  follow  each  particle  in  our 
simulation  for  At  seconds  according  to  the  equation  of  motion  for  a  single 
ion,  Eq.  (A-19). 

3.  Next  we  randomly  choose  a  collision  partner  for  each  ion. 

4.  Using  Monte  Carlo  method  we  then  determine  the  velocity  increments  for 
all  colliding  pairs  as  described  in  section  2.3.2.  The  new  velocities  are  fed 
back  into  the  Runge-Kutta  solver. 

5.  After  each  collision,  we  store  the  value  of  the  scattering  angle  0  for  each 
particle.  We  assume  that  whenever  sin2  ©  >  1  the  particle  has  under¬ 
gone  one  ion-ion  Coulomb  collision.  Here  the  summation  is  over  successive 
collisions  for  a  given  particle. 

2.3.2  Momentum  exchange  during  a  collision 

We  treat  Coulomb  collisions  between  ions  as  a  small  angle  binary  elastic  scat¬ 
tering  events  [27],  Such  collisions  preserve  energy  and  momentum. 

The  relative  velocity  vector  u(ux,uy,uz)  for  a  colliding  pair  is: 


U  =  Va  -  Vfc, 


(A-6) 
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where  va  and  v;,  are  the  velocities  of  two  colliding  ions.  The  post-collision 
relative  velocity  is: 

uf  =  \V  +  Au,  (A-7) 

where  u'  is  the  relative  velocity  right  before  the  collision  and  Au  is  the  change 
in  the  relative  velocity  due  to  the  collision.  The  change  in  the  relative  velocity 
due  to  a  scattering  event  could  be  derived  from  conservation  principles  [27], 

Aux  =  (ux/iix)u2sin©cos$  —  (uy/uj.)usin©sin$ 

-  ux(l— cos©),  (A-8) 

Aliy  =  (uy/u±)uzsmQ  cos4?  +  (ux/ux)usin©sin€> 

—  Uy(l  —  cos©),  (A-9) 

Auz  =  — ux  sin  ©  cos  $  —  nx(l  —  cos  ©),  (A-10) 

where  u2  =  u2+u2  and  u2  =  u2 +u2.  Here  perpendicular  and  parallel  directions 
are  defined  relative  to  the  magnetic  field  (i-axis).  When  ux  =  0  we  have, 


Aux 

=  u  sin  ©  cos  $, 

(A-ll) 

Ally 

=  usm©sin$, 

(A-12) 

Auz 

=  —  u(l  —  cos©). 

(A- 13) 

Angle  $  is  chosen  homogeneously  randomly  from  0  to  27T.  Angle  ©  is  chosen 
according  to: 


sin©  = 


1  —  cos  ©  = 


25 

1  +  <52  ’ 
2<52 

1  +  52’ 


(A-14) 

(A-15) 


where  5  —  tan(©/2)  is  a  random  number  chosen  with  the  Gaussian  distribution 
centered  around  zero  and  having  the  following  variance  ( 5 2): 


(. 5 2)  =  At 


q4n\ 

7reo77l2llull3 ' 


(A- 16) 


where  q  and  m  are  the  charge  and  the  mass  of  the  ion,  n  is  the  particle  number 
density,  A  is  the  Coulomb  logarithm,  eo  is  the  permittivity  of  free  space,  ||u|| 
is  relative  speed  of  two  colliding  ions,  and  At  is  the  time  step  [27].  This  small 
angle  restriction  allows  us  to  interpret  as  the  ion-ion  collision  frequency  vu  to  the 
binary  collision  frequency  iq,  and  the  formalism  implicitly  accounts  for  electron 
shielding.  We  did  not  however  track  the  electron  dynamics  as  electrons  are 
not  expected  to  be  effected  by  low  frequency  of  the  wave.  We  assume  that  the 
electron  temperature  is  not  affected  by  the  low  frequency  waves  of  the  problem, 

t O  ^  (jJq£  . 

The  post-collision  velocity  of  each  particle  is  found  simply  from, 

=  v„  +  (m/2)Au, 
vb  =  v6  +  (m/2)Au. 


(A-17) 

(A-18) 
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2.3.3  Moving  the  particles 

Starting  from  the  Lorentz  force  equation, 

F  =  mx  =  q( E  +  v  x  B),  (A-19) 

we  can  derive  equations  of  motion  for  a  single  particle  in  three  dimensions.  In 
our  analysis,  the  magnetic  field  is  constant  ,B  =  B z,  and  the  electric  field  arises 
from  the  propagating  electrostatic  waves,  as  shown  in  Fig.  2.1, 


y  +  E  ^2  sin(x  -  aijt), 

(A-20) 

X 

-x, 

(A-21) 

o, 

(A-22) 

so  that  x,jj,  and  z  are  the  second  derivatives  with  respect  to  time  t,  and  the 
other  variables  are  the  same  as  those  appearing  in  equations  (A-2)  and  (A- 
3).  Equations  (2.3.3)  could  be  solved  numerically  using  4th  order  Runge-Kutta 
method. 


2.4  Simulation 

The  above  model  is  used  to  simulated  the  case  of  BEW  acceleration  (e  = 
10,  iq  =  24.3,  i>2  —  25.3  m  —  1)  and  compared  to  those  of  SEW  acceler¬ 
ation  under  similar  condition  (e  =10,  u  =  24.3).  To  visualize  the  numerical 
results  we  use  Poincare  cross-sections  ( p  vs.  0  phase  diagrams)  [12].  We  will  also 
investigate  how  collisions  influence  the  energy  evolution  of  the  entire  system. 

2.4.1  Phase  diagrams 

Figure  2.3  follows  collisionless  evolution  of  1000  particles  in  the  plot  p  vs.  6. 
Initially  we  distribute  all  particles  homogeneously  over  region  of  the  phase  space 
p  <  20.  The  stochastic  heating  is  observed  whenever  ions  reach  the  stochastic 
zone  (p  >  20),  Fig2.2b.  Particles  with  initial  conditions  lying  outside  prohibited 
zone  are  accelerated  as  could  be  seen  from  that  figure.  The  points  corresponding 
to  unaccelerated  particles  define  a  mount-like  structure,  seen  in  the  last  two 
panels  of  Fig.  2.3,  which  corresponds  to  prohibited  zone  shown  in  figure  2.2b. 
Figure  2.4  shows  the  BEW  ion  acceleration  case  where  the  evolution  of  phase 
space  points  is  qualitatively  different  then  the  collisionless  case  illustrated  in 
2.3.  Even  ions  originally  in  the  forbidden  acceleration  zone  are  accelerated. 

2.4.2  Energy  Evolution 

Now  that  we  showed  that  collisions  enhance  ion  heating,  we  will  analyze  the 
energy  evolution  of  the  entire  system.  In  this  section  we  compare  the  cases 
of  BEW  and  SEW  ion  acceleration  by  beating  electrostatic  waves  and  a  single 
electrostatic  wave. 
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Figure  2.3:  £  =  10,  v\  =  24.3,  U2  =  25.3.  Collisionless  case.  The  ’’hump” 
corresponds  to  the  particles  not  accelerated  by  the  beating  waves  in  accordance 
with  Eq.  (2.2). 


Figure  2.5  shows  perpendicular  component  of  the  energy  for  BEW  as  well 
as  SEW  cases.  Initially  the  energy  increases  exponentially.  This  corresponds  to 
ions  funnelling  to  the  stochastic  region.  The  process  is  analogous  to  phase  space 
diffusion  -  thus  its  exponential  nature.  As  more  particles  find  their  way  into  the 
stochastic  region,  the  exponential  increase  is  followed  by  the  equilibration  of 
the  energy  (by  stochastic  motion  and  collisions).  Because  particles  fill  up  the 
stochastic  region  randomly,  the  statistical  average  of  the  energy  stays  constant. 

Figure  2.5  demonstrates  that  collisions  significantly  increase  both  the  ion 
heating  rate  and  the  final  average  energy. 

In  order  to  help  us  further  illustrate  these  effects  better  we  have  ran  the 
simulation  with  the  collision  frequency  restrained  to  a  constant  value  irrespec¬ 
tive  of  the  temperature.  Figure  2.6  shows  the  results  for  four  such  cases  and  a 
simulation  with  self-consistent  collision  frequency  calculation,  for  comparison. 
The  figure  indicates  that  there  is  an  optimum  collision  frequency  at  which  the 
heating  rate  and  efficiency  are  at  the  maximum.  When  the  collision  frequency 
is  increased  from  0  to  v\,  =  0.04wc,  the  ions  diffuse  faster  into  the  stochastic 
region  and  the  heating  efficiency  improves.  However,  as  the  collision  frequency 
is  further  increased  (yt>  =  4.14cjc),  the  heating  rate  as  well  as  the  efficiency  start 
dropping  because  collisions  increasingly  disturb  ion  motion  too  much.  In  a  real 
plasma  the  ion  collision  frequency  scales  as  ~  T~^2  (for  T i  <  Te).  As  the  elec¬ 
trostatic  waves  deposit  their  energy  into  the  plasma,  the  ions  collide  less  often 
and  the  collision  frequency  drops  to  its  optimum  value.  However,  if  the  collision 
frequency  decreases  even  further,  the  heating  efficiency  drops,  driving  down  the 
temperature  and  stabilizing  the  collision  frequency  back  to  its  optimum  value. 
Therefore,  we  can  conclude  that  the  collision  frequency  always  changes  to  ac- 
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Figure  2.4:  e  =  10,  v\  =  24.3,  =  25.3.  Particles  are  allowed  to  collide  with 

each  other.  We  take  the  initial  collision  frequency  required  for  step  1  of  our 
algorithm  to  be  ~  106,  which  corresponds  to  ne  ~  1012  cm-3  and  Te  =  300  K. 
Unlike  Fig.  2.3,  beating  electrostatic  waves  accelerate  all  ions. 


commodate  the  maximum  possible  heating  rate.  This  becomes  more  evident 
by  comparing  the  clamped  value  of  Vb  simulations  (dashed  lines)  to  the  solid 
line  which  was  obtained  by  running  the  simulation  and  allowing  the  collision 
frequency  to  change  self-consistently  with  the  ion  temperature  evolution. 

2.5  Conclusions 

Numerical  simulations  of  the  nonlinear  interaction  of  magnetized  ions  with  beat¬ 
ing  electrostatic  waves  (BEW)  were  carried  out.  The  resulting  particle  heating 
was  compared  to  that  obtained  from  simulations  with  interaction  with  SEW 
under  the  same  conditions.  The  higher  heating  rate  and  temperature  attained 
in  BEW  acceleration  can  be  explained  through  the  following  fundamental  de¬ 
scription.  In  the  SEW  interaction  thermalizing  collisions  are  the  only  means  for 
particles  below  a  certain  threshold  of  energy,  associated  with  the  resonant  con¬ 
dition,  to  reach  the  region  of  phase  space  were  stochastic  and  vigorous  heating 
takes  place.  The  non-resonant  character  of  BEW  acceleration  allows  a  signif¬ 
icant  fraction  of  the  ion  distribution  function  to  be  accelerated  and  reach  the 
stochastic  region.  This  acceleration  is  augmented  by  collisions.  In  addition  to 
the  thermalizing  role  of  collisions,  this  simulation-supported  phenomenological 
picture  points  to  the  promise  of  using  BEW  as  a  new  and  efficient  method  for 
accelerating  magnetized  ions  in  a  real  plasma. 
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Time  (seconds'coc) 

Figure  2.5:  Perpendicular  energy  evolution  for  1000  particles  interacting  with 
beating  waves,  e  =  10,  iq  =  24.3,  v2  =  25. 3. For  comparison  we  also  show  the 
energy  evolution  for  the  single  wave-ion  interaction,  e  =  10,  v  =  24.3. 


Time  (seconds*toc) 

Figure  2.6:  Perpendicular  energy  evolution  for  1000  particles  interacting  with 
beating  waves.  Dashed  lines  correspond  to  the  clamped  values  of  collision 
frequency  and  solid  curve  represents  the  self-consistent  simulation,  e  =  10, 
vi  =  24.3,  v2  =  25.3. 


Chapter  3 

Beating  Waves  Experiment: 
Excitation  of  Ion  Cyclotron 
Waves  and  Demonstration 
of  Ion  Acceleration 

3.1  Introduction 

Inductive  rf  plasma  sources,  and  in  particular  helicon  plasma  discharges,  are  of 
interest  in  propulsion  research  because  of  their  high  plasma  production  efficiency 
and  controllability  [28,  29].  Unfortunately,  these  rf  plasma  sources  produce  cold 
(Ti  0.1  eV)  ions.  To  achieve  high  Isp  it  is  thus  necessary  to  heat  these  ions 
considerably. 

Various  types  of  electrodeless  plasma  heating  (ion  acceleration)  provide  an 
efficient  way  to  increase  ion  temperature.  Methods  such  as  the  Ion  Cyclotron 
Range  Heating  (ICRH),  Lower  Hybrid  (LH)  wave  heating,  and  the  current  drive 
have  been  suggested  for  ion  heating  in  fusion  devices  [30,  31,  32].  The  ICRH 
scheme  is  also  employed  in  the  VASIMR  experimental  rocket  concept  [33,  34,  35]. 
In  addition,  ion  heating  by  various  electromagnetic  and  electrostatic  instabilities 
has  been  observed  in  the  Earth  ionosphere  [36,  37,  38,  39,  40]. 

In  this  chapter  we  investigate  excitation  of  Electrostatic  Ion  Cyclotron  (EIC) 
waves  that  propagate  transversely  to  the  external  magnetic  field.  Our  previous 
theoretical  and  numerical  studies  have  shown  that  two  beating  electrostatic 
waves,  obeying  specific  criteria  [10,  25,  41,  42],  may  energize  ions  very  effi¬ 
ciently.  Such  ion  energization  mechanism  can  prove  to  be  useful  in  propulsion 
applications. 

There  is  a  lack  of  detailed  measurements  of  EIC  wave  properties  and  prop¬ 
agation  in  rf-sust.ained  plasmas.  It  is  with  this  goal  in  mind  that  we  investigate 
the  excitation  and  propagation  of  EIC  waves  across  a  magnetized  rf-sustained 
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plasma  column,  and  subsequent  ion  energization  by  these  waves.  In  the  above 
sections  we  have  not  dealt  with  ion  energization  but  focused  on  the  excitation 
and  propagation  of  EIC  waves.  The  energization  is  demonstrated  in  the  last 
section. 

This  chapter  is  organized  as  follows.  In  section  3.2  we  review  previous  ex¬ 
perimental  work  on  electrostatic  wave  launching.  Then  in  section  3.3  and  3.4 
we  describe  our  experimental  apparatus  and  the  diagnostics  used  to  detect  and 
study  the  waves.  In  section  3.5  we  analyze  the  electrostatic  waves  launched  in 
our  apparatus.  We  summarize  our  results  regarding  wave  excitation  in  section 
3.7.  We  finally  conclude  in  the  last  section  with  the  experimental  demonstration 
of  the  existence  of  the  new  ion  acceleration  mechanism. 


3.2  Review  of  previous  work 

Various  experiments  reported  on  vigorous  ion  energization  by  EIC  waves  in  mag¬ 
netized  plasmas.  However,  most  of  these  studies  relied  on  exciting  EIC  waves 
through  some  internal  plasma  instability.  Two  typical  experimental  configura¬ 
tions  are  shown  in  Fig.  3.1. 

In  these  experiments  an  electrostatic  wave  of  frequency  ~  wCi  is  excited  ei¬ 
ther  by  drawing  electron  current  along  the  magnetic  field  to  a  positively  biased 
small  electrode,  as  shown  in  Fig.  3.1a,  or  by  creating  an  electric  field  perpen¬ 
dicular  to  the  external  magnetic  field,  as  shown  in  Fig.  3.1b.  In  both  cases 
significant  ion  energization  was  observed  once  the  wave  was  excited  [44,  45]. 
Unfortunately  in  experiments  like  these,  it  is  impossible  to  separate  the  cause 
from  the  effect  -  the  ion  energization  from  the  wave  generation  mechanism.  In 
addition  the  wave  frequency  cannot  be  controlled.  Thus  to  investigate  ion  ener¬ 
gization  properly  one  needs  to  design  an  experiment  where  the  waves  are  excited 
by  some  externally  controlled  antenna. 

Four  such  antenna  configurations  are  shown  in  Fig.  3.2.  Hooke  and  Barn- 
abei  [46]  used  capacitively  coupled  plates,  shown  in  Fig.  3.2a,  and  Stenzel  and 
Gekelman  [47]  employed  a  set  of  wires  strung  along  the  magnetic  field,  as  shown 
in  Fig.  3.2b,  to  launch  waves  close  to  the  LH  resonance.  Schmitt  launched  Pure 
Ion  Bernstein  Waves  (PIBW)  with  a  single  wire  at  the  center  of  the  plasma 
column  [48],  Fig.  3.2c,  while  Schmitt  and  Krumm  launched  Bernstein  waves 
through  a  mode  conversion  mechanism  with  a  wire  coil  wrapped  around  the 
plasma  column  [49],  Fig.  3. 2d. 

Goree  et  al.  [50]  and  Skiff  et  al.  [6]  used  electrostatic  plate  antenna  to  launch 
a  single  electrostatic  wave  above  the  ion  cyclotron  frequency  transversely  to  the 
magnetic  field.  Significant  stochastic  ion  energization  was  reported  in  the  latter 
experiment.  Since  our  theoretical  and  numerical  investigations  have  focused  on 
similar  frequency  range  we  have  adopted  this  antenna  design  for  our  experiment. 
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b) 

IS3  CSD  CS]  ESI  CSD 


ISD  CSJ  ISD  CS3  iS3 

Figure  3.1:  Typical  setup  for  studying  Electrostatic  Ion  Cyclotron  (EIC)  wave, 
a)  Used  by  Motley  and  D’Angelo  [43]  to  excite  current  driven  EIC.  The  picture 
is  taken  from  Ref.  [44],  b)  Used  by  Koepke  et  at.  to  excite  inhomogeneous 
energy  density  driven  EIC.  The  picture  is  taken  from  Ref.  [45]. 
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Figure  3.2:  Various  antenna  designs  for  launching  waves  into  a  plasma,  a) 
Capacitively  coupled  plates,  and  b)  a  set  of  wires  strung  along  the  magnetic 
field  launch  waves  close  to  the  LH  resonance,  c)  A  single  wire  at  the  center  of 
the  plasma  column,  and  d)  a  coil  wrapped  around  the  plasma  column  launch 
Bernstein  waves. 

3.3  Experimental  setup 

3.3.1  Vacuum  Chamber 

A  schematic  of  the  Beating  Wave  experimental  apparatus  (BWX)  is  shown  in 
Fig.  3.3.  It  consists  of  two  pyrex  cylinders  placed  inside  a  0.1  Tesla  magnet.  The 
axial  magnetic  field  along  the  centerline  is  shown  in  Fig.  3.4.  The  two  curves 
correspond  to  the  ion  cyclotron  frequency  of  10  kHz  and  30  kHz  in  the  test 
section  of  the  vacuum  chamber.  The  small  cylinder  is  6  cm  in  diameter  (ID) 
and  37  cm  in  length  while  the  large  cylinder  is  20  cm  in  diameter  (ID)  and  46  cm 
in  length.  The  backplate  of  the  small  cylinder  is  made  from  molybdenum  and 
is  electrically  floating  to  minimize  sputtering.  The  two  cylinders  are  connected 
by  an  electrically  floating  aluminum  plate  with  a  6  cm  concentric  hole  at  the 
center  to  allow  free  flow  of  gas  between  the  cylinders.  A  uniform  fill  pressure  of 
1  to  30  mTorr  is  maintained  by  a  gas  feed  (Ar  or  He)  at  the  aluminum  endplate 
of  the  large  cylinder  and  by  a  150  1/s  turbo  pump  with  a  conductance  controller 
backed  up  by  a  roughing  pump.  The  system  is  capable  of  maintaining  a  base 
pressure  of  2  •  10“6  Torr. 

Once  the  plasma  discharge  is  ignited  in  the  small  cylinder,  the  plasma  prop¬ 
agates  along  the  magnetic  field  lines,  which  are  parallel  to  the  axis  of  the  cylin¬ 
ders,  into  the  large  chamber  where  the  wave-launching  and  plasma-energization 
experiments  are  conducted. 


Pyrex  Cylinder: 
OD  =  22  cm 
ID  =  20  cm 
L=46cm 

Magnet 
Coils 

Figure  3.3:  The  drawing  of  the  BWX  experimental  apparatus. 
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Figure  3.4:  The  axial  magnetic  field  (Bz)  along  the  centerline  of  the  magnet. 
The  ion  cyclotron  frequency  is  relatively  constant  within  the  test  portion  of  the 
large  chamber  where  the  electrostatic  waves  are  launched. 
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Figure  3.5:  Plasma  density  and  electron  temperature  for  various  rf  power  to 
the  saddle  antenna.  As  rf  power  is  raised  above  500  W  the  high  plasma  density 
helicon  discharge  is  observed.  B=782  Gauss,  P=lmTorr. 


3.3.2  Plasma  Source 

A  Boswell  saddle  type  antenna  used  to  create  the  plasma  discharge  is  placed 
around  the  small  cylinder.  The  antenna  is  made  of  0.25”  copper  tubbing  to 
allow  water  cooling.  An  inductive  discharge  is  produced  by  supplying  rf  power 
to  the  antenna  from  an  ENI  13.56  MHz  1.2  kW  power  supply  through  a  tuner. 
The  tuner  consists  of  an  L  network  made  of  two  Jennings  1000  pF  3  kV  variable 
vacuum  capacitors.  The  tuner  is  placed  as  close  to  the  antenna  as  possible  to 
maximize  coupling. 

An  inductive  discharge  is  easily  obtained  with  only  a  few  watts  of  forward 
rf  power  to  the  antenna  and  only  a  few  percent  of  rf  power  reflected.  A  helicon 
discharge  with  high  plasma  density  (1013  cm-3)  can  be  produced  as  rf  power  to 
the  antenna  is  raised  above  500  W  by  properly  adjusting  the  pressure  and  mag¬ 
netic  field.  The  inductively  coupled  discharge  looks  homogeneous  and  occupies 
the  entire  cross-section  of  the  small  cylinder.  When  the  helicon  discharge  is  ob¬ 
tained  with  argon,  a  bright  blue  column  is  observed  at  the  centerline.  Measured 
radial  distributions  of  plasma  density  and  electron  temperature  for  both  types 
of  the  discharges  are  shown  in  Fig.  3.5.  Figure  3.6  shows  the  range  of  some 
natural  plasma  frequencies  that  were  calculated  with  measured  Te,  ne,  and  B. 

3.3.3  Beating  Waves  Antenna 

In  attempt  to  launch  electrostatic  waves  into  the  plasma  column  we  have  tried  a 
couple  of  antenna  configurations.  An  external  antenna  consisting  of  two  spools 
of  wire  in  the  Helmholtz  coil  configuration,  placed  on  opposite  sides  of  the  out- 
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Figure  3.6:  Typical  plasma  parameters  for  the  BWX  experiment  calculated  with 
measured  Te,  ne. 


side  wall  of  the  vessel,  with  their  axis  perpendicular  to  the  magnetic  field  were 
initially  used  but  did  not  excite  electrostatic  waves.  The  hope  was  that  electro¬ 
magnetic  waves  from  this  external  antenna  could  be  converted  into  electrostatic 
waves  by  the  plasma  [51,  52].  We  were  able  to  launch  an  electrostatic  wave 
with  an  antenna  placed  inside  the  vacuum  chamber,  and  consisting  of  two  flat 
metal  plates.  This  type  of  antenna  has  been  used  to  launch  a  single  electrostatic 
wave  above  the  ion  cyclotron  frequency  by  Goree  et  al.  in  the  toroidal  ACT-1 
device  [50,  53],  and  by  Skiff  et  al.  in  a  linear  device  [6].  The  antenna  is  made  of 
two  1  cm  x  6  cm  molybdenum  plates  placed  3.14  cm  apart  along  the  magnetic 
field.  The  plates  are  oriented  such  that  the  long  side  and  the  surface  normal 
are  perpendicular  to  the  magnetic  field,  as  shown  in  Fig.  3.3.  During  the  wave¬ 
launching  experiments  reported  below,  the  plates  were  driven  by  a  Wavetek  180 
signal  generator  with  a  sinusoidal  signal  through  a  Tektronix  AM  501  modular 
op-amp  either  in  or  out  of  phase,  as  shown  in  Fig.  3.7. 


3.4  Diagnostics 

Two  types  of  Langmuir  probes  were  employed  to  measure  the  steady-state 
plasma  properties  and  the  wave  propagation.  Plasma  density  and  electron  tem¬ 
perature  were  determined  using  a  radio  frequency  compensated  Langmuir  probe 
with  0.5  mm  graphite  tip  [54].  The  rf  compensation  was  achieved  with  four 
miniature  inductors  placed  in  series  and  close  to  the  probe  tip.  The  inductors 
where  chosen  to  filter  out  the  fundamental  at  13.56  MHz  and  the  second  har¬ 
monic  at  27.12  MHz  of  the  helicon  antenna  signal.  Uncompensated  Langmuir 
probes  with  either  graphite  or  T-shaped  tips  made  of  10  mil  tungsten  wire  were 
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Figure  3.7:  The  circuit  diagram  for  launching  a  single  electrostatic  wave  into 
the  plasma  column.  The  two  antenna  plates  can  be  driven  in  or  out  of  phase. 


Figure  3.8:  To  launch  a  pair  of  beating  waves  the  impedance  between  the  two 
signal  sources  and  the  amplifier  has  to  be  matched  through  a  power  combiner. 
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used  to  detect  the  wave  and  investigate  its  dispersion. 

Instead  of  the  Laframboise  analysis  [55]  of  the  Langmuir  probe  I-V  charac¬ 
teristic  we  used  an  empirical  floating  potential  method  specifically  developed  by 
Chen  et  al.  [56]  to  determine  Te  and  ne  with  rf-compensated  Langmuir  probes 
in  helicon-sustained  plasmas  similar  to  ours. 


3.5  Wave  launching 

The  circuit  used  to  launch  a  single  electrostatic  wave  into  the  plasma  column 
is  shown  in  Fig.  3.7.  The  sinusoidal  signal  (10-300  kHz)  from  a  Wavetek  180 
signal  generator  is  amplified  with  gain  of  10  by  a  Tektronix  AM  501  modular 
operational  amplifier  (Rout;  ~  150  SI,  ±40  V,  50  mA  Max.).  To  drive  the  plates 
in  phase  the  signal  is  sent  to  both  plates.  A  coupling  capacitor  (Cc  =  2pF)  can 
be  used,  as  indicated  in  Fig.  3.7,  to  allow  both  plates  to  float  with  respect  to 
the  plasma  potential.  To  drive  the  plates  out  of  phase  the  antenna  is  connected 
through  a  1:1  transformer. 

The  launching  of  a  pair  of  beating  waves  into  the  plasma  is  accomplished  by 
combining  the  output  from  two  signal  sources.  This  is  achieved  by  mixing  the 
low-power  signals  coming  out  from  two  Wavetek  signal  generators  and  sending 
the  combined  output  to  an  amplifier,  as  demonstrated  in  Fig.  3.8.  The  output 
impedance  of  the  signal  generators  and  the  input  impedance  of  the  amplifier  is 
50  fl.  To  make  sure  that  the  signal  of  one  of  the  signal  generators  is  not  distorted 
by  the  other,  we  needed  to  match  the  impedance  of  the  entire  circuit.  A  power 
combiner  consisting  of  a  100  fl  resistor  and  a  tapped  inductor  in  parallel  provide 
good  matching  between  the  signal  generators  and  the  amplifier. 

To  ensure  that  the  waves  are  launched  into  the  plasma  efficiently  we  also 
needed  to  maintain  a  good  antenna-plasma  coupling.  Different  amount  of  power 
is  delivered  to  the  antenna  plates  depending  on  the  impedance  mismatch  be¬ 
tween  the  driving  antenna  circuit  and  the  plasma.  Plasma  impedance  varies 
with  the  plasma  density,  i.e.  with  rf  power  delivered  to  the  helicon  antenna, 
and  the  fill-up  pressure.  For  a  given  pressure  the  power  to  the  antenna  plates 
and  the  efficiency  of  the  electrostatic  wave  launching  can  be  correlated  to  the 
helicon  antenna  power. 

In  order  to  investigate  the  wave  launching  conditions  we  conducted  the  fol¬ 
lowing  experiment.  The  signal  generator  and  the  op-amp  in  Fig.  3.7  were  set  to 
produce  a  30- Volt  sinusoidal  signal  (uj  =  70  kHz)  with  no  plasma.  A  magnetized 
argon  plasma  with  ne  =  109  — 1013  cm-3  and  fcj  =  30  kHz  was  produced  by  the 
helicon  antenna  with  the  chamber  back  filled  to  1  mTorr.  Measuring  the  plates 
voltage  and  current  simultaneously  we  were  able  to  determine  the  impedance 
and  the  power  delivered  to  the  antenna  plates  as  a  function  of  rf  power  deliv¬ 
ered  to  the  helicon  antenna.  A  Langmuir  probe  inserted  into  the  plasma  column 
measured  the  amplitude  of  the  launched  waves. 

We  found  that  for  the  in-phase  configuration  the  maximum  power  to  the 
antenna  plates  is  delivered  when  rf  power  is  175-200  W,  as  shown  in  Fig.  3.9a  and 
3.9b.  The  amplitude  of  the  Langmuir  probe  signal  is  a  complicated  function  of  rf 
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Figure  3.9:  The  amplitude  of  the  launched  wave  as  well  as  the  impedance  and 
the  power  delivered  to  the  antenna  plates  are  shown  as  functions  of  the  rf  power 
to  the  helicon  antenna.  The  real  and  imaginary  components  of  the  antenna 
impedance  are  designated  by  Zr  and  Z*  respectively.  The  fundamental  as  well 
as  the  second  harmonic  of  the  launched  waves  are  measured  with  a  Langmuir 
probe,  a)  Both  plates  are  driven  in  phase  (without  the  coupling  capacitor),  b) 
Both  plates  are  driven  in  phase  (with  a  coupling  capacitor  connected  in  series 
with  the  antenna,  Fig.  3.7).  c)  The  plates  are  driven  out  of  phase. 
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power,  and  it  does  not  seem  to  track  the  power  delivered  to  the  antenna  plates. 
There  are  two  maxima  for  the  detected  wave  amplitude.  One  occurs  at  200  — 
250  W  of  rf  power,  when  the  real  component  of  the  measured  impedance  ( ZT ) 
approaches  150  fi,  the  output  impedance  of  the  op-amp,  and  the  imaginary  part 
( Zi )  is  small.  At  that  condition  the  entire  circuit  is  matched.  The  second,  and 
higher  maximum  occurs  at  350—400  W  of  rf  power,  and  probably  corresponds  to 
the  increase  in  plasma  density  due  to  the  transition  from  the  inductive  discharge 
mode  to  the  helicon  mode.  The  real  part  of  the  measured  impedance  at  that 
point  is  50  fl.  A  more  efficient  wave  launching  can  thus  be  achieved  by  choosing 
an  amplifier  with  an  impedance  close  to  that  value. 

It  is  interesting  to  note  that  the  Langmuir  probe  detected  a  plasma  wave 
at  the  second  harmonic  aside  the  fundamental  frequency  driven  by  the  circuit. 
As  shown  in  Fig.  3.9  the  second  harmonic  can  be  a  significant  fraction  of  the 
fundamental. 

A  similar  experiment  was  repeated  with  the  antenna  plates  driven  out  of 
phase.  The  results  are  shown  in  Fig.  3.9c.  The  maximum  power  delivered  to 
the  antenna  in  this  case  is  about  75%  of  the  power  in  the  case  described  above. 
Also,  the  wave  amplitude  detected  by  the  Langmuir  probe  is  significantly  lower 
than  the  amplitude  of  the  wave  launched  by  the  plates  driven  in  phase.  Possible 
reasons  as  to  why  the  out-of-phase  configuration  does  not  couple  well  into  the 
plasma  will  be  given  in  the  next  section.  However  as  in  the  previous  case,  the 
maximum  in  the  wave  amplitude  corresponds  to  the  matched  circuit  condition, 
Zi  ZT  ~  150  n,  and  not  to  the  maximum  power  delivered  to  the  antenna 
plates.  The  imaginary  component  of  the  impedance  for  the  out-of-phase  con¬ 
figuration  is  negative,  indicating  capacitive  coupling  between  the  plates.  In  the 
in-phase  configuration  the  imaginary  component  is  positive,  indicating  inductive 
coupling. 

The  in-phase  driven  antenna  couples  better  to  the  plasma  thus  resulting 
in  a  wave  with  a  higher  amplitude  than  the  out-of-phase  driven  antenna.  We 
therefore  undertook  further  study  of  the  in-phase  configuration.  The  goal  of 
the  study  was  to  determine  whether  the  launched  wave  is  electrostatic  in  nature 
and  whether  it  propagates  transversely  to  the  magnetic  field. 


3.6  Dispersion  relation  measurement 


The  electrostatic  dispersion  relation  derived  from  the  kinetic  theory  for  an  in¬ 
finite  slab  of  collisionless,  homogeneous,  isotropic,  and  non-drifting  plasma  can 
be  written  as  [50,  30], 

D  =  k\Kxx  +  k\Kzz  =  0,  (A-l) 

where  k±  and  fc||  are  the  perpendicular  and  the  parallel  components  of  the 
wavenumber,  and  Kxx  and  Kzz  are  two  components  of  the  dielectric  tensor. 
For  a  single  ion  species  plasma  these  components  can  be  expressed  as, 

Kxx  —  1  + 
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where  In  is  the  modified  Bessel  function,  A j  =  k\Tj/m,jU)2j,  Z(//)  is  the  plasma 
dispersion  function  [57],  and  (i  =  (w  —  niVj)(mj /2T,)1/2 /fcjj,  with  j  =  e,i.  In 
Eqs.  (A-2)Eq:Kzz  we  assumed  that  w  <8C  f le  and  Ae  4C  1. 

To  launch  an  electrostatic  ion  cyclotron  (EIC)  wave  above  the  ion  cyclotron 
frequency  the  following  conditions  must  also  be  satisfied  [50], 


bJpx  >  ujcz  ,  T{  5^  Te, 

(2Ti/mi)1/2  <  w/*y  <  (2 Te/me)1/2. 

In  our  experiment  we  have  Te  =  3  eV,  Tt  =  0.1  eV,  B  =  261  Gauss,  ui  = 
30  -  180  kHz,  and  Ay  «  46  cm,  therefore  the  above  inequalities  are  satisfied. 
Here  we  assumed  that  the  parallel  wavelength  is  determined  by  the  extent  of  the 
large  glass  cylinder  where  the  waves  are  launched.  By  measuring  the  parallel 
component  of  the  wavenumber  with  two  axially  separated  Langmuir  probes 
we  have  confirmed  that  the  parallel  wavenumber  is  ~  46  cm.  The  numerical 
solution  of  Eq.  (A-l)  for  the  plasma  parameters  in  our  experiment  is  represented 
in  Fig.  3.10. 

Theoretical  curves  can  be  subdivided  into  two  types  according  to  their  slopes. 
One  curve  extending  from  k±  =  0  to  2  cm-1  has  a  positive  slope.  This  curve 
represents  the  forward  branch  of  the  EIC  dispersion  relation  since  its  group 
velocity  ( dw/dk )  is  in  the  same  direction  as  its  phase  velocity  (oj/k).  On  the 
other  hand,  curves  extending  from  k±  =  1  to  6  cm-1  have  negative  slopes.  The 
curves  correspond  to  cyclotron  harmonics  of  the  backward  branch  of  the  EIC 
dispersion  relation,  since  their  group  velocity  has  the  opposite  sign  of  the  phase 
velocity.  The  pass-band  of  each  of  these  backward  branches  is  narrow,  and 
therefore  we  do  not  expect  to  see  any  of  the  backward  mode  waves  propagating 
in  our  experiment. 

The  slope  of  the  forward  branch  of  the  EIC  dispersion  is  sensitive  to  the 
electron  temperature,  while  the  backward  branch  is  sensitive  to  the  ion  tem¬ 
perature.  Thus  comparing  the  experimentally  obtained  dispersion  relation  to 
the  theoretical  expression  provides  a  good  check  of  the  species  temperatures. 
While  the  ion  temperature  was  not  measured,  the  electron  temperature  inferred 
by  the  slope  of  the  experimental  dispersion  relation,  Te  =  2.75  eV,  is  in  good 
agreement  with  the  values  measured  independently  with  the  rf-compensated 
Langmuir  probe,  shown  in  Fig.  3.5. 

As  was  mentioned  above,  the  theoretical  dispersion  relation  given  by  Eq. 
(A-l)  was  derived  for  an  idealized  plasma.  Some  plasma  parameters  of  our 
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Figure  3.10:  The  measured  dispersion  relation  at  fcj=10  kHz  agrees  well  with  the 
theoretical  dispersion  for  the  forward  branch  of  the  Electrostatic  Ion  Cyclotron 
wave  (EIC).  Curves  extending  to  large  values  of  k±  represent  the  cyclotron 
harmonics  of  the  EIC  backward  branch. 

experiment  are  given  by  Fig.  3.6.  It  can  be  seen  that  the  plasma  is  collisional 
(z/ej7  Ua  >  fCi).  Also,  the  density  profiles  shown  in  Fig.  3.5  indicate  that  the 
plasma  is  not  homogeneous,  (a/n)dn/dx  >  1,  specifically  at  high  rf  power. 
Another  point  of  concern  is  the  effect  of  the  plasma  boundaries.  At  low  values 
of  k±  (long  wavelength)  the  effect  of  the  boundaries  should  be  significant  and 
the  experimental  data  might  diverge  from  the  theoretical  expression.  With 
these  considerations  in  mind,  Fig.  3.10  shows  a  surprisingly  good  agreement 
between  the  experiments  and  the  theory  for  a  wide  range  of  frequencies.  The 
experimental  data  were  obtained  in  the  following  manner. 

A  single  wave  was  launched  by  the  antenna  plates  inserted  at  the  edge  of  the 
plasma  column  and  driven  in  the  in-phase  configuration.  We  measured  the  wave 
dispersion  by  a  system  of  three  uncompensated  Langmuir  probes.  The  probes 
were  placed  orthogonally  to  each  other  so  that  simultaneous  measurements  of  k± 
and  fc||  could  be  performed  by  measuring  the  signal  delay  between  any  two  probes 
[58].  Measuring  the  probe  signal  with  the  plasma  discharge  off  we  determined 
that  the  ac-coupling  signal  was  two  orders  of  magnitude  below  the  wave  signal 
with  the  plasma  turned  on. 

The  parallel  component  of  the  wavenumber  measured  by  the  two  probes 
placed  along  the  magnetic  field  indicated  that  the  parallel  wavelength  of  the 
wave  is  46  cm  -  the  length  of  the  large  glass  cylinder.  This  measurement  was 
confirmed  by  the  wave  observations  at  various  axial  positions  along  the  chamber. 
The  perpendicular  wavenumber,  shown  in  Fig.  3.10,  varies  in  accordance  with 
the  EIC  dispersion.  The  dispersion  measurements  indicate  that  the  wave  is 
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propagating  at  the  angle  of  82°  to  86°  with  respect  to  the  magnetic  field.  Both 
the  fundamental  and  the  second  harmonic  components  of  the  signal  were  used 
to  construct  the  experimental  dispersion  relation  shown  in  Fig.  3.10. 

Driving  the  plates  of  the  antenna  out  of  phase,  we  expected  the  parallel 
wavelength  to  be  6.28  cm  -  twice  the  spacing  between  the  plates.  This  however 
was  not  the  case.  The  measured  parallel  wavelength  was  on  the  order  of  the 
chamber  length.  The  measurements  of  the  perpendicular  wavelength  were  not 
consistent  with  the  EIC  dispersion  relation.  It  is  likely  that  the  wave  was  affected 
by  the  ion  Landau  damping  since  w/k^  ~  vthi- 

3.7  Conclusions  regarding  ES  Wave  Excitation 

We  reported.on  the  excitation  and  propagation  of  an  Electrostatic  Ion  Cyclotron 
(EIC)  wave  launched  by  a  two-plate  antenna  into  a  magnetized  argon  plasma. 
We  described  in  detail  the  circuity  necessary  to  drive  the  antenna  in  the  in-  and 
out-of-phase  configuration. 

We  found  that  the  antenna  driven  in  phase  couples  better  to  the  plasma  and 
excites  higher  amplitude  waves  than  the  plates  driven  out  of  phase. 

In  addition,  we  determined  that  the  wave  amplitude  can  be  optimized  by 
carefully  choosing  the  right  plasma  parameters  and  building  a  matched  driving 
circuit. 

The  antenna  driven  in  phase  excites  a  wave  with  Ay  =  46  cm  -  the  length  of 
the  portion  of  the  vacuum  chamber  where  the  waves  are  launched.  The  perpen¬ 
dicular  component  of  the  wavenumber  varies  according  to  the  EIC  dispersion 
relation  and  is  in  good  agreement  with  the  forward  branch  of  the  theoretical 
dispersion  relation  despite  the  simplifying  assumption  of  a  collisionless,  mag¬ 
netized,  homogeneous,  isotropic,  and  infinite  plasma  slab.  We  found  that  EIC 
waves  launched  at  frequencies  between  ojci  and  I0wci  propagate  with  little  damp¬ 
ing  at  an  angle  between  82°  and  86°  with  respect  to  the  magnetic  field. 

Multiple  EIC  waves  required  for  our  ongoing  beating  waves  studies  [10,  25, 
26,  21,  42]  can  also  be  easily  launched  by  the  plates  antenna  with  two  or  more 
signal  generators  by  means  of  a  power  combiner. 

3.8  Demonstration  of  Ion  Acceleration 

After  the  success  in  launching  the  waves,  we  proceeded  to  measure  the  effects 
of  the  beating  waves  on  ion  acceleration  to  provide  the  first  test  for  the  ex¬ 
istence  of  the  new  ion  acceleration  mechanism.  The  goal  was  to  measure  the 
ion’s  perpendicular  enregy  (equivalently  velocity  or  temperature)  and  observe 
if  ion  acceleration  (or  equivalently  heating)  occurs  when  a  pair  of  beating  ion 
cyclotron  waves  are  launched.  In  order  to  do  this  measurement  we  developed, 
using  funds  from  a  recent  DURIP  grant,  a  dedicated  state-of-the-art  Laser  In¬ 
duced  Fluorescence  (LIF)  diagnostic  system  (shown  in  Fig.  (3.11))  that  allows 
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Figure  3.11:  Photograph  and  schematic  of  the  new  state-of-the-art  LIF  system, 
recently  developed  under  DURIP  funding  and  used  for  our  AFOSR-funded  beat¬ 
ing  wave  study. 


accurate  measurement  of  the  ions’  energy  distribution  function.  Our  LIF  system 
is  described  in  a  recent  publication  ??. 

The  measurements  led  to  the  recently  obtained  results  plotted  in  Fig.  (3.12) 
which  shows  as  much  as  40%  enhancement  to  the  ion  energy  due  to  beating 
wave  acceleration.  This  represents  the  first  experimental  evidence  ever 
of  this  new  acceleration  mechanism. 


Figure  3.12:  Measured  ion  energies  (blue  data  points)  under  the  effect  of  a  pair 
of  beating  ion  cyclotron  waves  compared  to  the  benchmark  case  with  no  waves 
(red  horizontal  line  at  .1  eV)  showing  that  the  ions  were  accelerated  by  as  much 
as  40%.  (RF  Power  in  helicon  source  =  250  W,  Wave  Antenna  voltage  =  60  V, 
ion  cyclotron  frequency  =  10  kHz,  wave  frequency  =  10-100  kHz.) 
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Impact  Statement 


The  revised  budget  represents  about  two-thirds  the  original  budget.  Therefore,  the  scope 
of  the  work  will  be  adjusted  as  follows. 

The  original  plan  was  to  explore  the  use  of  three  different  layered  materials:  bismuth 
telluride,  bismuth  selenide,  and  molybdenum  diselenide.  However,  in  view  of  the 
reduced  budget,  we  will  concentrate  on  the  first  and  last  as  the  first  two  are  quite  similar. 
I  do  not  think  this  will  seriously  affect  the  scientific  understanding  of  nanocomposites  we 
expect  from  this  study. 


